Spin Related Effects in Transport Properties of "Open" Quantum Dots 
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We study the interaction corrections to the transport coefficients in open quantum dots (i.e. dots 
connected to leads of large conductance G 3> e 2 /irh), via a quantum kinetic equation approach. 
The effects of all the channels of the universal (in the Random Matrix Theory sense) interaction 
Hamiltonian are accounted for at one loop approximation. For the electrical conductance we find 
that even though the magnitude of the triplet channel interaction is smaller than the charging energy, 
the differential conductance at small bias is greatly affected by this interaction. Furthermore, the 
application of a magnetic field can significantly change the conductance due to the Zeeman splitting, 
producing finite bias anomalies. For the thermal conductance we find that the Wiedemann-Franz 
law is violated by the interaction corrections, and we investigated the effect of magnetic field on 
the Lorentz ratio for contacts of finite reflection. The charge and triplet channel corrections to 
the electrical and thermal conductance vanish for reflectionless contacts. In the latter case the 
temperature and magnetic field dependence of the conductance is determined by the Maki- Thompson 
correction in the Cooper channel. 

PACS numbers: 73.23.-b,73.21.La,73.63.Kv,71.70.Gm 



I. INTRODUCTION 

Transport properties of quantum dots are strongly af- 
fected by electron-electron interactions, the most stud- 
ied example being the so-called Coulomb blockade 
phenomenon!* 2 ^ (see Ref. for a review) in dots con- 
nected to external electrodes by contacts of low conduc- 
tance, G <C e 2 /irh = G q : if the temperature T is smaller 
than the charging energy E c associated with the addi- 
tion of one electron to the dot, transport through the 
dot is exponentially suppressed (except at charge degen- 
eracy points). As the contacts' conductance increases, 
quantum fluctuations of the dot's charge will eventually 
lift the Coulomb blockade; even so, interactions in such 
"open" (G> G q ) or weakly blockaded dots do influence 
transport phenomena^*^ In this paper we consider the 
interaction corrections to transport coefficients of first or- 
der in 1/g = Gq/G, similar to the well-known corrections 
in the bulki& 

It is known that transport properties of dots with a 
large number of electrons can be described within the 
framew ork o f Random Matrix Theory (RMT) (see e.g. 
Refs. I"il l9ll()l for reviews). We consider a system sepa- 
rated into clean leads without interaction, a closed dot 
and the coupling between the dot and the leads as shown 
in Fig. n All interference and interaction effects are as- 
sociated with the closed dot Hamiltonian Hrj. In the 
RMT approach, the orbital dynamics of non-interacting 
spinless electrons in the closed dot is described by the 
M — > oo limit of an M x M Hamiltonian matrix with 
random entries H^ nn belonging to the Gaussian ensem- 
ble: 




x=0 



FIG. 1: Quantum dot connected to two leads, and the corre- 
sponding Hamiltonians. 



gh is a dimensionless parameters quantifying the orbital 
effect of the weak magnetic field. This can be estimated 
as: 
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where Si is the one-particle mean level spacing, (• • •)„„ 
denotes averaging over the random matrix ensemble, and 



where Et is the Thouless energy! (above which scale the 
RMT description is not applicable), <& is the total mag- 
netic flux through the dot, and <E>o — hc/e is the flux 
quantum. The extreme values, g = and g = AM, cor- 
respond to the orthogonal and unitary ensembles respec- 
tively, while intermediate values describe the crossover 
between those. 

In the absence of spin-orbit interactions (we neglect 
spin-orbit interactions in this paper; the interplay of 
exc hang e and spin-orbit interactions was studied in 
Ref.[llD, the spin of the electrons can be accounted for by 
defining the following non-interacting Hamiltonian oper- 
ator (in the second quantized notation): 

H° = 1&£On, (1-3) 

where ("0m) is a t wo component spinor operator, 
whose components create (annihilate) electrons in the 
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m-th orbital state, with specific spin projections. For 
interacting electrons, the closed dot Hamiltonian is given 
by 



Hd — H° + H int , 



(1.4) 



where the dominant part of the interaction Hamiltonian 
has the universal form4*i^ 



H int = E C N 2 + jj 2 + J C T^T. 



Here 



A I 



n=l 
1 M 

n=l 
1 M 



(1.5) 

(1.6a) 
(1.6b) 
(1.6c) 



are respectively the number of electrons in the dot, the 
total spin of the dot electrons, and the pairing operator. 
Here, a and a v are Pauli matrices acting in the spinor 
space. 

The last term in the right hand side of Eq. (|1.5|) 
describes pairing between electrons (interaction in the 
Cooper channel) and for J c < it drives the dot towards 
the superconducting state; studies of such superconduct- 
ing grains have been reviewed e.g. in Ref. ll3L Here we 
will assume that either J c > and hence no supercon- 
ducting transition at any temperature, or J c < but 
T — T c > T c , so that the dot is in the normal phase, and 
furthermore, the superconducting fluctuations are small. 
The second term (triplet channel) gives the dependence 
of the dot's energy on the total spin in the dot - the ef- 
fects of this term on the tunneling density of states and 
on the spin susceptibility' 1 ' 4 and on the peak spacing 15 
in the Coulomb blockade regime have been recently con- 
sidered. Finally the first term (singlet channel) describes 
the charging energy and it is responsible for the Coulomb 
blockade. In the weakly blockaded regime, only this term 
has been considered previously in the literature* 5 *^ as 
its contribution is expected to be the dominant one (for 
the repulsive Coulomb interaction). However this term 
is not affected by an external magnetic field, whereas the 
two remaining terms are. Our goal is to calculate the 
interaction corrections to the transport coefficients with 
the full universal Hamiltonian Ho taken into account 
and to examine the dependence of these corrections on 
the applied magnetic field. In particular, we consider 
the non-linear conductance for voltage-biased dots and 
their (linear response) thermal conductance. To evaluate 
these transport coefficients we construct, starting from 
the RMT description, a quantum kinetic equation anal- 
ogous to the one developed for the description of disor- 
dered metals 



The remainder of the paper is organized as follows: in 
the next section we summarize our results for the conduc- 
tance of metallic quantum dots. In Sec. 1 1 1 1 1 we present 
the derivation of the O-dimensional Usadel equation in 
RMT. This equation is the starting point for the deriva- 
tion of the kinetic equation as outlined in Sec. II VI In 
Sec. we give the explicit calculation of the transport 
coefficients. 



II. SUMMARY OF THE RESULTS 

Here we present our results for the interaction correc- 
tions to the differential electrical conductance and the 
linear thermal conductance of quantum dots in the pres- 
ence of a magnetic field. These results are derived in 
Sec. El 



A. Electrical Conductance 

The total ensemble averaged differential conductance 
G is 



AG, 



where 



Go 



gLgR 



irk g L + g R 



(2.1) 



(2.2) 



is the classical conductance and AG is the interaction 
correction. We do not include the weak localization cor- 
rection for non-interacting electrons which can be found 
e.g. in Ref. 0. We also do not study the contribu- 
tion of AG to the mesoscopic fluctuations of the con- 
ductance, which are smaller than AG by the additional 
factor G q /Go, (For the charge channel this effect was 
studied in Ref. |?J). 

In Eq. (|2.2|l g£ and gR are the dimensionless conduc- 
tances of the left and right contacts respectively: 



T 
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T 
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(2.3) 



where T n is the transmission coefficient of the n-th chan- 
nel. 

The interaction correction AG has distinct contribu- 
tions from each term in the interaction Hamiltonian (|1.5f) : 



AG = AG C + AG S + AG. 



Cooper ■ 



(2.4) 



We postpone the discussion of the Cooper channel cor- 
rection, AGcooper, until the end of this section. For the 
charge and triplet channel contributions we have found 



AG C 



e 2 h L g| + h R g 2 L n 
irh (gL + gfl) : 



L E(T Q ,AE C ,T;V). 



(2.5) 
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and 

nh (g L + g K ) 3 

x ^ H(r + im^,j s ,r ; y). 

rn=0,±l 

The form factors h^and h^j are given by 

neL(R) 



(2.6) 



(2.7) 



Factors of this form were first obtained in Ref. [la and 
reproduced in Refs. This structure was originally 

missed in the formalism of Ref. 5 but was recovered in 
Ref.0 

In Eqs. (|2.5|) and l|2.6|) the dimensionless function S is 
defined by 



S(I\£,T; V) 



Re- 



where 



T-ieV 
2nT 



9(z) = V (0) («) + ^ w («) 
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(2.9) 



and (z) is the i-th derivative of the digamma function. 
The form of this function agrees with the result derived 
in Refs. 00 for the charge channel. We denote with V 
and T the bias voltage and the temperature respectively, 
Tq/H is the escape rate: 



h 6 Nch 
r o = - = — > 

n=l 



(2.10) 



and E* z is the Zeeman energy renormalized by the ex- 
change interaction: 



E* z = 



9l^bB 

1 + J B /5i ~ 1 + J s /S x 



(2.11) 



Here B is the magnetic field, is the Lande g-factor, 
and \ib is the Bohr magneton. 

In Fig. H we plot AG = AG C + AG S in units of 
G q = e 2 /irh, for different values of the triplet channel 
interaction constant J s in the absence of magnetic field. 
We have taken Nl = Nr and T n = 1/2 for all channels, 
so that the prefactors in Eqs. (|2.5|l and H2.6|) are equal to 
1/8. The interaction constant takes values ranging from 
J s = to J s = — 0.75<5i, which apply to most metals and 
quantum dots (the values of J s /8\ are reviewed e.g. in 
Ref. Irif) . The charging energy E c on the other hand is 
much larger than the mean level spacing. A theoretical 
estimate for two dimensional dots 4 yields E c /5± ~ r s kpL, 
where r s is the gas parameter, fej? is the Fermi wave- 
length, and L is the lateral dimension of the dot. This 
is large because of the large factor kpL 3> 1. As shown 
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FIG. 2: Interaction correction to the dimensionless differential 
conductance for varying strength of the triplet interaction in 
the absence of Zeeman splitting. Here T = 0.1<5i, To = 5<5i 
and E c — 100<5i where <5i is the mean level spacing, and we 
have taken Nl = Nr and T„ = 1/2 for all channels. The 
value of J 3 is shown above each graph. For clarity, graphs are 
shifted upwards by 0.05 at each step. 
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FIG. 3: The magneto-conductance of the starred curve of 
Fig.[2]vs. bias voltage , for different values of Zeeman splitting 
energy (shown above each graph). The graphs are shifted 
upwards by 0.03 at each step. 



in Fig. the charge channel correction generates a dip 
in the conductance at zero bias with a characteristic 
width of order (l + 4E c /Si)To, analogous to the zero bias 
anomaly in bulk systems. This width can also be esti- 
mated as h/r c , where r c < r is the classical recharging 
time governing the charge dynamics. On the other hand, 
the attractive triplet interaction produces a peak with a 
much smaller width of order (1 + J s /5i)Tq = H/t s , cor- 
responding to the slow spin dynamics: t s > r. Therefore 
we have a competition between the two corrections at 
zero bias as seen in Fig. [21 

In Fig. we plot the magneto-conductance, G(Ez) — 
G{Ez — 0), for different values of the Zeeman energy 
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FIG. 4: The magneto-conductance of the starred curve of 
Pig- vs. Zeeman splitting energy, at various bias voltages 
(shown above each graph). The graphs are not shifted here. 



(and with the same assumptions about the form factors 
made for Fig. [21 . In the presence of a magnetic field, the 
triplet channel contribution decomposes into three terms 
due to the Zeeman splitting. These terms produce peaks 
at eV ~ ±E% and 0. As the value of the Zeeman energy 
increases so does the width of the displaced peaks (for 
m = ±1); in the limit Ez 3> To this width is given by 
the bare Zeeman energy Ez- In Fig. 0] we plot AG vs. 
Ez for different bias voltages. 

In the case of refiectionless contacts, (T„ = 1 for 
all n), the corrections in Eqs. I|2.5|) and (|2.6(l vanish 
due to the vanishing of the form factors h^ and h^ of 
Eq. I|2.7|) : a non-zero contribution to G is given by the 
zero-dimensional analog of the Maki-Thompson 19 - 20-21 
correction to the conductivity: 



AG, 



Cooper 



AG 



MT 
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gLgR 



E* ) TTh (gL +gfi) 

x Yl S a g^ a0 (T„E* z ,T;V) 

a.fi=L,R 



(2-12) 



where T ll = T^/j and Tlr. = Trl are dimensionless 
functions given below, and T*/h is the escape rate mod- 
ified by the effect of the the magnetic field on the orbital 
motion of the electrons [see the discussion after Eq. 



Here 



r = n _ Sx 

1 * — — c, 
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'N ch 
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n=l 



max{T, eV,r»}' 



(2-13) 



(2-14) 



and T c is defined in Eq. 1|4.1U2[1 . In the attractive case, T c 
is the critical temperature of the superconducting tran- 
sition. In this case we only consider the normal state at 
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FIG. 5: Correction to dimensionless differential conductance 
for refiectionless contacts and symmetric leads at zero tem- 
perature and r„ = 55i. Here we have plotted £ 2 AGmt, so 
that the actual magnitude of the correction is further sup- 
pressed by 1/e 2 . The values of the effective Zeeman splitting 
E z are shown below each graph, and the graphs are shifted 
upwards by 0.25 at each step. 



Et S> T ^> T c , and hence |e| ^> 1. In the repulsive case, 
T c ^> Et, and since by assumption Ex is larger than 
all the other energy scales in the problem, e>1, This 
means that the correction in Eq. (|2.12|l is logarithmically 
suppressed for normal dots and is in general much smaller 
than the correction due to charge and triplet channels. 
For this reason we report here the Cooper channel cor- 
rection only for refiectionless contacts, where it is the 
only non-vanishing correction. The results for the gen- 
eral transmission coefficient can be found in Sec. IV CI 

The functions T aj g cannot be calculated in a compact 
form and we therefore consider approximate expressions 
valid in certain regions of parameters. In the low tem- 
perature, low voltage regime (i.e. T, eV <C I\) we have: 



T 



a/3 



[l + (E* z /T^] 2 



A a p 



eVy f_(T_ 

(2.15) 

where Alr. = Arl = 1, and All = Arm = 1/4. The 
T 2 and V 2 dependences are due to the inelastic nature 
of the processes at the origin of the Maki- Thompson 
correction 

At finite bias and low temperature (T <C eV, r*), we 
have: 



T 



LR 



LL 



- E h 

2 ^-f 

m=±l 



eV eV-mE% 



F, 



(2.16) 




arctan ■ 



eV - mE% 



m=±l 



where the full expression for function h is given in 
Eq. (|5.64|) . This function is used in Fig. [S] to plot 
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e 2 AGmt (in units of e 2 /ttK) for symmetric leads at dif- 
ferent values of the Zeeman energy. An approximate ex- 
pression for h is: 



h 



E 



arctan 



arctan 



2eV + 2(m - n)£| 



(2.17) 



eV - 



r, 



— arctan 



r* 



This formula qualitatively renders the shape of the cor- 
rection, although it overestimates it in the central plateau 
[cf. Fig[S]; on the other hand, Eq. i|2.17|) is a good ap- 
proximation when \V\ > E* z and E* z > 1/t*, i.e. it gives 
a quantitative description of both the steps at V ~ ±J5J 
and the asymptotic regions at large bias for sufficiently 
large Zeeman energy. At small bias (V <C 1/t*), the 
correction is better represented by Eqs. 12.15fl . 

Finally in the high temperature regime (T r») we 
have: 



LR 



• coth 



Ez_ 
2T 



T 



LL 



where 



E 



C2 



m— ± 

mE"; 



E% 



meV 



2T 



2T 



c n (x) 



d n 
dx n 



(x coth x) 



(2.18) 



(2.19) 



B. Thermal Conductance 

For the thermal conductance we calculate the linear 
response in singlet and triplet channels only, and we do 
not report the contribution of Cooper channel, as it is 
smaller by the factor 1/e 2 [cf. Eq. (|2.14|l ]. 

The thermal conductance is a combination of two parts 



KWF + Ak. 



(2.20) 



T L -T R 

The first term k,wf respects the Wiedemann-Franz law: 



n W F = f^G(V = 0) = ^ [Go + AG(V = 0)] , 

(2.21) 

with AG given by Eqs. I|2.5fl - H2.6|) . whereas the correc- 
tion Ak violates this law: 



Ak = Ak c + Ak s 



Ak c 



g|h L +g|h fl irT 
(gi + gi?) 3 9fi 
( 2ttT 
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(2.22a) 



(2.22b) 
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FIG. 6: The relative change in the Lorentz number vs. 
temperature for different values of Zeeman splitting energy 
(shown above each graph). Here, E c = 100<5i and J s = 
— 0.45<5i, To = 5<5i, and we have taken Nl = Nr = 1 and 
T n = 1/2 for all channels. The graphs are shifted upwards by 
0.03 at each step. 



Ak s = 



g 2 R h L + g 2 L h R TlT 
(gL + gi?) 3 9ft 



Re £ 



r 



2ttT 



9i 



,(l + lf) r - 



9i 



m=0,±l 

In Eqs. 12.22fl we use the notation 

r m = T + imE* z , 

and 



9i (x) 



(1) 



6 



3 

x ' 



(2.22c) 
/2ttT n 

v t; 

(2.23) 
(2.24) 



where ij)^ 1 ' is the derivative of the digamma function. The 
deviation from the Wiedemann-Franz law can be quanti- 
fied by defining a generalized Lorentz number according 
to 



L 



TG(V = 0)' 



(2.25) 



such that it would yield the usual L — it 2 /3e 2 in the 
absence of interactions. In Fig. [5] we plot the relative 
change in the Lorentz number, L/Lq — 1, as a function 
of temperature for different values of Zeeman splitting 
energy (and with the same assumptions about the form 
factors made for Fig. EI> • 

We note that the interaction corrections to thermal 
conductance Eqs. (|2.22b|) - (|2.22c|l vanish for reflection- 
less contacts, and so the Wiedemann-Franz law is satis- 
fied in this case. That this law is not violated by inelastic 
processes in this case can be explained by the following 
qualitative picture. Consider an electron with energy e 
and charge e entering the dot from the right lead with an 
influx proportional to the number of channels N R in this 
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lead (T n = 1 for all channels). The electron will subse- 
quently leave the dot through the right or the left lead 
with probability N R /N c h and Nl/N^ respectively, thus 
contributing an outflux in the right lead proportional to 
Nfi/Nch- Therefore, the total ingoing electric current in 
this lead will be 



I R cx eN R 



c N 2 R _ N L N R 
6 N ch 6 N l +Nr 



(2.26) 



i.e. the electric conductance is given by the classical con- 
ductance Eq. H2.2fl . If in the process the electron loses an 
energy oj to the collective excitations this will still hold 
within the linear spectrum approximation as the bosonic 
excitations are neutral and do not contribute to electric 
current. This is apparently not the case for the energy 
current as the electron going out of the right lead now 
has energy e — oj. However, the collective excitation with 
energy u> will eventually decay into an electron hole pair 
and these will enter the leads with the same probabilites 
as before, and due to energy conservation will carry the 
lost energy back to the right lead with the same rate, i.e. 



1% oc eN R 



(e-co) 



N 2 R 



N 2 R 



N, 



ch 



N 



ch 



e NlNr . (2.27) 
N L + N R 1 ; 



We see that the relation between the energy current and 
the electric one Eq. I|2.26[l is not altered and we recover 
the Wiedemann-Franz law. It is crucial in this reasoning 
that the contacts be reflectionless, since only in this case 
the transmission coefficients will not be renormalized by 
interaction and the T„'s remain equal to one independent 
of energy. This is not the case for T n ^ 1 and the renor- 
malized transmission coefficients depend on energy, thus 
violating Eq. (|2.27|l and therefore the Wiedemann-Franz 
law. 



III. DERIVATION OF THE OD USADEL 
EQUATION 

In this section we derive the zero dimensional Us- 
adel equation^ that describes electrons in open quan- 
tum dots, taking into account the electron-electron inter- 
action but neglecting weak localization and mesoscopic 
fluctuation effects. Our treatment here is based on the 
Keldysh technique^ for non-equilibrium systems, and is 
analogous to Ref. ^3 with appropriate modifications for 
the 0-d case. First we start with the description of cou- 
pling to the leads in subsection IIII Al In the rest of the 
subsections we derive equations for Green functions av- 
eraged over the appropriate random matrix ensemble. 
To overcome the difficulties involving ensemble averag- 
ing in the presence of the quartic interactions of the 
universal Hamiltonian Eq. (|1.5I) . we employ a Hubbard- 
Stratonovitch transformation for each channel of interac- 
tion, rendering the fermion Hamiltonian quadratic. This 
is done in subsections IIII Bl to IIII D| 



A. Description of the Leads 

The open dot [see Fig. ^ is described by the Hamil- 
tonian^ 



H = H T 



H L + H 



LD, 



(3.1) 



where H D = H° + H int is the (RMT) Hamiltonian for 
the interacting electrons in the dot [see Eqs. i|l.lfl - l|1.5|) ]. 
Hl is the Hamiltonian for the free electrons in the leads 
and Hld describes the coupling between the dot and the 
leads. The electron spectrum in the leads near the Fermi 
surface can be linearized: 



f dk 

H L =v F Y,J 2^i(*)V>a(*), 



(3.2) 



where vp is the Fermi velocity, and a labels different 
channels in the leads: 1 < a < Nl for the left, and 
Nl + 1 < a < N c h for the right lead channels, and N c h = 
Nl + N R . Here, the field ip a is understood to be a two 
component spinor and we suppress the spin indices unless 
stated otherwise. 

The leads-dot coupling Hamiltonian is 



Hld= ]T e-^iWn^Uk^n + h.C.), 



where the coupling constants are defined as 



tr 



M5i 

o Van i 



(3.3) 



(3.4) 



for n ~ a < N c h and zero otherwise, and the exponential 
at 77 — > + , is used to regularize the coupling at large \k\. 
Here v = \/{2itvf) is the density of states per spin at the 
Fermi level. We can always write the matrix W in the 
above diagonal form by choosing the appropriate basis 
for the random matrix. 

We introduce the exact Green functions of the elec- 
trons in the dot (G) and in the leads (F). As usual they 
are 2x2 matrices in the Keldysh and spin spaces 



n ( 17 nm 

(h;t 2 ) QtL(h;t 2 ) 

r , FKJtxM-MM) F^Jt 1 ,k 1 ;t 2 ,k 2 ) 
P 1 FL{tiM\t 2 ,k 2 ) Ffjti,ki;t 2 ,k 2 ) 



(3.5) 



K 



where 1 < n,m < M (M — > oo being the size of the 
random matrix), and a, [3 label different channels in each 
lead. The entries of the matrix in Eq. I|3.5|) are given by 

F R = -i6i 2 (ip a (k u ti)V»J(fca, h) + ^(k 2 ,t 2 )ip a {k u h)), 
F A = i9 2 x{ip a (ki,ti)ipl(k 2 ,t 2 ) +ipl(k 2 ,t 2 )ip a (k 1 ,t 1 )), 



= -i(i/j a (h,ti)ipl(h,t 2 ) - t/j0(k 2) t 2 )ip a (ki,ti)), 



F z 



0. 



(3.6) 



The expression for Q is obtained from Eq. I|3.6|) by re- 
placing a, (3 with n, m and removing the k variables. 
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Here 9 12 = 6{t\ — t 2 ), where 0(t) is the step function, 
the fermionic spinor operators are in the Heisenberg rep- 
resentaton and the averaging is performed over a non- 
equilibrium state of the system. 

For Hld = 0, one easily finds the Green functions for 
the electrons in the leads F = F^ ' 



F<-V = 2ir5 aP 6{k 1 - k 2 )h 

R 1 



de 

2V 



-^-^F^ick); 



4 0) (e,fc) 



e — vpk + iO ' 
1 

e — vpk — iO ' 



K 



F<?\e,k) =-irid(e-v F k)fK(e), 



(3.7) 



where l s is the unit matrix in the spin space, and 
—2 < f K (e) < 2 is an arbitrary function related to the 
occupation number, n(e), of the state with energy e by 
71(e) = j — t/(c)- For the equilibrium Fermi distribution, 
/(e) is given by 



can write the expectation value of the outgoing current 
in each lead in terms of the Keldysh Green function 



Ia{t) 



-evp 



dkdq sin rjiq r 



Tt s F K {t,k + q;t,q), 



(3.12) 

where a = L, R for the left or right lead, a is summed over 
different channels in the corresponding lead, and Tr s is 
trace in the spin space. After substituting Eq. (|3.7|) into 
Eq. H3.9fl . to obtain F , and integrating over momenta 
we obtain 



In 



2tt 



A ) 

II I! J 



Here the first term in the parentheses has the meaning 
of the current from the dot to the leads, and the other 
terms express the current from the leads to the dot. Note 
however, that the terms taken separately do not have any 
physical significance. 



/f (e) = 2tanh 



2T„ 



(3.8) 



In the presence of coupling with the dot, we can use the 
standard diagrammatic technique [see Fig. 0d] to write 
the following expression for F 

F af} (k 1 ,k 2 ) = F^(k 1 ,k 2 ) + f |«e-^M^ 



x F^>(h,k[)W a , n g nm Wl p J(%{k' 2 , k 2 ), (3.9) 

where summation over repeated indices is implied. We 
generally regard the Green functions as operators in the 
time as well as in the Keldysh and spin spaces, so that 
their products are understood as operator (matrix) mul- 
tiplication in these spaces, and therefore we omit time 
arguments in particular. 

Equation 13.9(1 enables us to analyze the transport 
properties of the system, e.g. the electric current. The 
linearized spectrum of the leads Hamiltonian Eq. (|3.2|) 
suggests a similar linearization for the current operator 
in the channel a (we choose the direction such that out- 
going current is positive) 



Ia{x) = —ipl l (x)d x i/} a (x) + h.c. 
Am 



(3.10) 



Substituting tp a ( x ) — f^l V'aCO sin(fc + kp)x, x < 
into Eq. I(3.1()|l . and keeping the terms not oscillating on 
the scale of the Fermi wavelength 2ir/kp, we find (for 

X = Til) 

aEa 

where ryx — > +0, 771 ^> r/, and integration is performed 
for |fc|,|g| <§C k F . Therefore, according to Eq. ((3.6|) we 



B. Charge Channel 

The charge channel of the universal interaction Hamil- 
tonian Eq. 1)1. 5J1 is given by 



(3.14) 



Flint — E C N 2 



where the number operator N was defined in Eq. 1(1. 6a() . 
This interaction can be decoupled via the Hubbard- 
Stratonovitch transformation by introducing a scalar 
fluctuating field <f> with the following matrix structure 
in the Keldysh space: 



(3.15) 



K 



This field modifies the self-energy of the electrons by con- 
tributing a term proportional to <f> [see Eq. 1(3.20(1 and 
Fig. 0b] . For the purpose of the one loop approxima- 
tion we will employ in the next section, these fields can 
be considered as Gaussian with the propagators 



+(*l)< 


Ufa))* 


= 5^(<l,*2) ) 


(3.16a) 


+ (*!> 


M*2)>* 


=\v*{ti,t 2 )> 


(3.16b) 


+ (*a)« 

-(«a)< 


M*i)>* 
M*0>* 


= l -vj{ti,t 2 ), 

=0. 


(3.16c) 
(3.16d) 



In the saddle point approximation, the propagators are 
solutions of the matrix (in Keldysh space) Dyson equa- 
tion 



= 2E, 



(i + n 
2>M 



(3.17) 



n„ 
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id. 



b) 



Qn 



= 1 + ( J 

4> w*> 



F'-°> " 
X-----X + (RMJ 



- + 



n m n n m. n m n 



The crossing diagrams in Fig. 0c are smaller than the 
non-crossing diagram by a factor of 1/M, and will be 
neglected. In addition, the second term in con- 
tributes to the non-crossing diagram only when m = n 
(see Fig. 0c), and is therefore small by a factor of 1/M. 
This means that the orthogonal and unitary ensembles 
are equivalent for charge and triplet channel (the Cooper 
channel case is considered in Sec. IIII Dl . Thus, in the 
M — > oo limit, we obtain the following Dyson equation 
for Q n 



F 



(e) 



I * 



-t 



9 W 
-x — «— x- 



FIG. 7: Dyson equation for the ensemble averaged dot Green 
function (a-c), and the expression for the lead Green function 
(d). The crossing diagrams in (c) are smaller than the non- 
crossing one by powers of 1/M and can be ignored. The 
disordered matrix elements in (e) are given by Eq. 



where 1 = lx5(ti — t 2 ), and 

n ( tl ,t 2 ) - n (t 2 , h)-^ • 



M 



nf(ii,i 2 ) = £ 



2i5<t>-(t 2 ) 



(3.18) 



Here we have introduced the Green function of the dot 
electrons as a functional of the field 4> hi its matrix form 
in Keldysh space: 

Qnmitum- {g^MMLfrM*))*' (3 - 19) 

such that its average over the fluctuating field 4> gives 
the usual expressions as in Eq. (|3.6[) . with the averaged 
Q z vanishing. We will suppress the argument <j> in the 
subsequent formulas. 

We construct perturbation theory for the RM aver- 
aged Green function (Q n m) RM = QnSnm, as an expansion 
in powers of the random matrix Hamiltonian Hq, and 
the Hubbard-Stratonovitch field <j> (see Fig. for the di- 
agrammatic form of the Dyson equation^* for Q n ). The 
first term in the self energy expansion (Fig. 0b) is due 
to the effect of the Hubbard-Stratonovitch field <j>, which 
only changes the phase of the Green function and does 
not lead to scattering between different orbitals. The 
second term is due to scattering to and from the leads. 
The third term (Fig. 0c) gives the contribution of the 
non-interacting closed dot Hamiltonian H° after averag- 
ing over the random matrix ensemble, where the ran- 
dom matrix correlation in panel (e) is given by l|l.l|l . 



. dGn c 



iwW na f a Wl n + M 



(3.20) 



2 M 



Here, as in subsection IIII Al the products are understood 
as operator (matrix) multiplication in time as well as in 
Keldysh and spin spaces. In particular for the unit oper- 
ator we have 1 = 8{t\ — t 2 )lK <8> l s , and <j> is understood 
as the operator <j>{t\)8{ti — t2)l s - We have also defined 



a = — [ ^e-^F^(k), (3.21) 
ttv J 2tt 



/< 



where Pi°\k) was defined in Eq. lpT7jl , Using Eq. p77Jl 
we integrate over k obtaining 



where 



fa{h-t 2 ) 



1 fS W 
-1 



(3.22) 



(3.23) 



K 



and the function /^(e) was defined in Eq. (|3~%|) . Substi- 
tuting Eq. (|3.4|) into Eq. (|3.20() and defining 



we obtain 

1 dg n 
M dh 



iM-Q n , 

7T 



1 



M „ 

E 9n 
M' 

n—1 



99n 



1, 



(3.24) 



(3.25) 



where we have adopted units such that Si/tt = 1: 
t 



f 



Si/tt 



Si 

7T 



As we are interested in the limit M — > oo of the RMT, 
we can neglect the first two terms on the left side of 
Eq. I3.25f> . obtaining 



\tn\ 2 Ugn + gg n = 1, (n < N ch ), 
gg n = i, (n > N ch ). 



(3.26) 
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Summing Eq. I|3.26[l over all n and neglecting terms of 
order N c h/M we obtain the following constraint for g: 

9-9 = 1. (3.27) 
Using Eqs. I|3.26[l - i|3.27[l we solve for g n 's, obtaining 



9n 



1 



g+\tn\ 2 fn 
(i + \tn\ 2 9fn) 



9(l + \tn\ 2 fn9Y 



(3.28) 



Thus, given g and the lead Green functions Eq. l|3.23|l . 
we can completely determine the g n 's using Eq. 1|3.28[1 . 
As for g, Eq. (|3.27|) in general only constrains g to a cer- 
tain manifold, and does not determine it further. The 
non-equilibrium evolution of g on this manifold is given 
by the terms of order 1/M in Eq. (|3.25(l which were ne- 
glected in reducing that equation to Eq. (|3.26(l . In order 
to seperate those 1/M terms, we subtract from Eq. (|3.25|) 
its transpose, obtaining 



1 d ~ 9n _l • 1 U - 1 



\tn\ 2 [fn \9n] + [ 9n ! 9\ , 



which after summation over all n gives 



E 



I L n \ 2 [fn ;g n 



(3.29) 



Here d/dt = d/dti + d/dt2 is the derivative with respect 
to the "center of mass" time, and [a;b] — ab — ba is the 
commutator. Using Eq. I|3.26[) . we can rewrite Eq. 
as 



dg 
dt 



■i[4>;g] 



E 

n=l 



(3.30) 



which we call the 0D Usadel equation^ It is easy to 
see that this equation is consistent with the constraint 
Eq. KTli . 

We are now prepared to take the M — * oo limit in 
the formula for the current Eq. Ij3.13|l . Using the defini- 
tion Eq. jSIllil and constraint Eq. (|3~2T)| in Eq. (|3~T3|l we 

obtain 



/ » = iE 1Y '(-i f »i 2 [-'»^J K )- 



?5> 



s [g ; 9n\ 



lA' 



(3.31) 



Thus the total outgoing current is proportional to the 
Keldysh part of the right hand side of the Usadel equa- 
tion, Eq. 1)3.30(1 . which states nothing but the conserva- 
tion of number of particles. 

We can rewrite the right hand side of Eq. (|3.3U|I in 
a form that is more convenient for later use. Using 



Eq. (|3.28() . we write this commutator in terms of g and 
fn only: 

[g\9n] = (i + M 2 /„.grMi + M'ffAr 1 
= \tn\ 2 [g;l 



9n ■ 9n 



1 + ^T n ({£;/„} -2 



(3.32) 

-l 



where {a; b} = ab + ba is the anti-commutator, and in 
obtaining the third line from the second, we have used 
Eq. !j3~2"8ll and Eq. iEO?) to write 

9n-g n = [(l + \t n \ 2 f n g) ■ (i + \tn\ 2 gfn)]~ l 

= "(l + l^| 2 ) 2 i + |tn| 2 ({5;/n}-2)] _1 . 

Here the transmission coefficient of the n-th channel is 
given by 

4|t n | 2 



(3.33) 



To obtain the polarization operators, we average 
Eq. (|3~T%|) over realizations of the RMS Using Eq. (533^ . 
we obtain 



nj(ti,ta) 



nf(ti,t 2 ) 



-_S(t 1 -t 2 ) 



254>+{t 2 ) 
STi- s (g K (t 1 ,t 1 )+g z (t 1 ,t 1 )) 



(3.34) 



28(1)- (ta) 

The appearance of the 8- function term here deserves com- 
ment. It is analogous to the usual ultraviolet anomaly. 
Namely, in the first line of Eq. (|3.34|) the expression 
g K (ti,ti) is understood as the limit t\ — > t[ of g K (ti, t[), 
where the latter is calculated in the M — > oo limit al- 
ready. On the other hand, Eq. (|3.18() implies the oppo- 
site order of limits, and the 5- function takes care of this 
discrepancy: 



-~5(t! 

IT 



t2) 



(3.35) 



-, M 



lim lim 



lim lim 



n=l 



2i5(j>+(t 2 



This equality can be checked by explicit calculation. 
In fact, the coefficient in front of the delta function is 
guarded by the requirement of gauge invariance, i.e. the 
total electron number response to the potential <j) u should 
vanish for u> 1/t. 

Equation (|3.30() together with the constraint l|3.27|l . 
the identity (|3.32|) . the propagators Eq. 1)3. 17|) and the 
polarization operators Eq. I|3.34[) . completely determine 
the kinetics of the quantum dot, given the distribution 
functions of the electrons in the leads. Moreover, the last 
equality in Eq. (|3.32|) gives a convenient starting point for 
the perturbation theory; indeed in the non-interacting 
case, g has the same structure as /„, Eq. (|3.23|l . and the 
square bracket in the last line of Eq. I|3.32|l reduces to 
the identity (and the commutator gives the relaxation of 
g towards /„). 
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C. Triplet Channel 

The triplet part of the interaction Eq. I|1.5() is given by 



(3.36) 



where the total spin operator S was defined in Eq. (|1.6b(l . 
This can similarly be decoupled using a vector Hubbard- 

Stratonovitch field h, with each component having the 
same Keldysh structure as in (|3.15() . The presence of 
this field leads to the replacement 4> —> <fi + h, 



h 



(3.37) 



in the self-energy Eq. (|3.20|l . Here i,j — x,y,z, and a 1 
are the Pauli matrices in spin space. The propagators 
for this field are defined as in Eq. (|3.16() by replacing 
the ^>'s with different components of the vector field h, 
such that now each Keldysh component of the propagator 
acquires a 3 x 3 tensor structure. In the saddle point 
approximation we have 



V h = 



(l + UhD h ) 



(3.38) 



where 1 = lif<5y<5(ii — t 2 ), and the polarization operator 
tensors are given similarly as in Eq. (|3.18|l by variational 
derivative of the Green functions with respect to com- 
ponents of the field h. Repeating the steps leading to 
Eq. I|3.34|l . we obtain 



[nfL(ti.ta) 



7T 



25hJ{t 
6Tr s (a*(g K +g z )(t 1 ,t 1 )) 
26hJ(t 2 ) 



(3.39) 



where i,j — x,y,z. Furthermore, in the presence of a 
magnetic field in the z-direction, we must add to the 
dot's Hamiltonian the Zeeman energy term 



H z = E z z-S, (E z = g L[ i B B), 



(3.40) 



where gi, is the Lande g-factor and fig is the Bohr mag- 
neton. This causes the ^-component of the field h to 
acquire a non-zero average at the saddle point, given by 



(hf) 



2 ' 



(3.41) 



where the other components of the tensor vanish due to 
the symmetry of spin rotations around the z-axis. Here, 
the factor 1/2 arises from the same factor accompany- 
ing the Pauli matrix in the definition of the operator S, 
Eq. (|1.6b|) . Separating this average from the fluctuating 



part, and redefining h to stand for the latter, we obtain 
the modified Usadel equation: 



dg_ 

dt 



-iE* z [a z ;g]+i[h;g]=J2 [ffS&J, (3-42 



where E* z is the renormalized Zeeman energy 



E* Z = E Z + 2(hi 



(3.43) 



Equation (ETUI together with Eqs. (|3~3S|) - (j3~3§|) . and 

the constraint 1)3.27(1 and identity Eq. (|3.32|1 which are 
still valid, completely determine the kinetics of the quan- 
tum dot under the effect of spin fluctuations. 



D. Cooper Channel 



Interaction in the Cooper channel is given by 



Hint = JcT^T 



(3.44) 



where the pairing operator T was defined in Eq. (|1.6c[l . 
This interaction can be decoupled using a complex 
Hubbard-Stratonovitch field A which also has the struc- 
ture l|3.15(l in Keldysh space. This pairing field gives rise 
to anomalous Green functions which can be taken into 
account by introducing the standard Gor'kov-Nambu 26 
(GN) spinor 



(3.45) 



where a y acts in the spin space, and we denote Pauli ma- 
trices in the GN space by t % . We also redefine the dot and 
lead Green functions Q nm and F a p to be 2 x 2 matrices 
in the GN, as well as in the Keldysh and spin spaces (we 
will use bold symbols to represent matrices in GN space) , 
such that after averaging over the Hubbard-Stratonovitch 
field different Keldysh components will have the same 
structure as in Eq. Q3.6f) . with ip and ip' replaced by t z ^ 
and ^ respectively, and the products understood as di- 
rect products in GN as well as in spin spaces. Because of 
the standard convention to include r z in the definition of 
GN Green functions, the time derivatives in the properly 
modified Eqs. I|3.20|l and (|3.30|l will be accompanied by 
this Pauli matrix [see l|3.53ll below]. We introduce the 
notation 



g = 




(3.46) 



to represent the GN components of the Green function 
g (related to £„'s as in Eq. $T2Q ). The upper left GN 
component of this matrix gives the Green function used 
in the singlet and triplet channels, and the off-diagonal 
components are the anomalous Green functions. Since 
the components of the GN spinor are not independent 
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= <j v ® t v ^), the components of the Green function 
g are not entirely independent either, and are related by 

g = Sg T S, (S = a%®a v ®t x ), (3.47) 

where af^ is the Pauli matrix acting in the Keldysh space, 
and the superscript T stands for transposition of the op- 
erator g in time, Keldysh, spin and GN spaces. The same 
relation is satisfied by /„'s so that their explicit form is 



f = I /n 

" ' fl rr 1 f T rr x , 

u a K Jn a K/ N 



(3.48) 



and the matrix /„'s are given by Eq. 1|3.23|) . We define 
the matrix A and the one-particle Hamiltonian H° to 
have the following forms in GN space: 



A 

-A* 



o 



Y 



H° 
a v H 0T a y 



(3.49) 



JV 



where we have included the Zceman splitting Hamilto- 
nian Eq. (|3.4U|) in H° . With these definitions the expres- 
sion for self-energy l|3.20|) takes the form (in units where 

h/it = i) 

S„ = k + ^--i-KvW n J a Wl n -iMg+i^T, (3.50) 
where we have defined 



E z = E* z a> 



LAY 



(3.51) 



and 



For repulsive interaction ( J c > 0) or for attractive in- 
teraction (J c < 0) at T > T c , the field A can be consid- 
ered as Gaussian with propagators 



(A + (1)A;(2)) A =-Z>£(l,2), (3.54) 

(A + (1)A*_ (2)^ = 1^(1,2), 

(A_(1)A;(2)) A = 1^(1,2), 
<A_(1)A*_(2)) A =0. 

In the saddle point approximation, the propagators sat- 
isfy the matrix (in Keldysh space) Dyson equation 



x>a = J c (i + n A r> A j , 



(3.55) 



where 1 = 1^-5(^1 — t 2 ), and the polarization operators 
are given by 



n2( f „« = ni( f2 , (l ;- iaT '-'^< f - f '» 



4 SA+(t 2 ) 



nf(ti,ta) = - A 



4 5A-(h 



(3.56) 



Equation 13.5311 together with Eqs. (|3.55|) - (|3.56|l . and 
the constraint l|3.27l) and identity Eq. (|3.32() (which are 
still valid with the additional GN structure understood), 
constitute a complete description of the kinetics of the 
quantum dot under the effect of pairing fluctuations. 



T = 




(3.52) 



The origin of the last term in Eq. (|3.5U|) is explained 
as follows: in the presence of a time reversal symme- 
try breaking mag netic field, H° and H 0T in Eq. (|3~4^|) 
are not equal. Therefore, to leading order in 1/M, the 
non-crossing diagram of Fig. 0c (being proportional to 
(H°gH°) RM ) receives contributions from the second term 
in the right hand side of Eq. i|l.l|) in the off-diagonal GN 
parts. 

Using Eq. (|3.5U|) and following steps similar to the ones 
leading to Eq. lEHUj) . we arrive at the properly modified 
Usadel equation for the GN Green functions 



dq dq , t . 

— + — — T + - 

dh dt 2 2 1 



^ Z l^- + ^ + -[E z ;g}-^[T;g} 



Sh, 
4 



n=l 



(3.53) 



We note that in order to take into account the interaction 
in all three channels together, we can simply add the 
fields = ® ljv and h = h ® t z to A in Eq. I|3.53|l . 



IV. DERIVATION OF THE KINETIC 
EQUATION 

Even though the Usadel equation l|3.53|l gives a com- 
plete description of the kinetics of the open quantum 
dot, its solution for a general form of the Hubbard- 
Stratonovitch fields is not tractable due to its non- 
linearity and non-locality in the time domain. The pur- 
pose of this section is to reduce the Usadel equation to 
the kinetic equation for all three channels of the interac- 
tion. The kinetic equation describes processes character- 
ized by the time At which are slow. In the particular case 
of the quantum dots, slowness means At 3> 1/w*, where 
lj* is the characteristic energy transfer in the electron- 
electron interaction, (the scale uo* is determined by the 
shape of the distribution function f K (e), e.g. for thermal 
equilibrium to* ~ T) . The relation between At and the 
relaxation scale determined by the kinetic equation itself 
(i.e. 1/r defined in Eq. (|2.10() ') may be arbitrary (see e.g. 
Sec. II B of Ref. [H for more detailed discussion). 

In the first loop approximation for the interactions, the 
different channels are decoupled and we will treat them 
seperately. 
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A. Charge Channel 

1. Gauge Transformation 

The condition that justifies the first loop approxima- 
tion in spite of the large charging energy (E c 3> S±), is the 
largeness of the leads' conductance which suppresses the 
effect of fluctuations of the Hubbard-Stratonovitch field 
for small frequency ui ~ 1/r ^> S\. On the other hand, 
if the transmitted frequency is larger than 1/t, the dot 
can be considered as closed. For the closed dot, the field 
<fi is coupled only with the total number of electrons, N, 
which commutes with the Hamiltonian. This interaction 
causes the motion of the one-particle energy levels with- 
out any redistribution of particles between them, which 
even though produces a large, singular effect in the one- 
particle Green function (known as Coulomb blockade), 
has nothing to do with relaxation inside the dot. To elim- 
inate the effect of the trivial motion of the levels from the 
very beginning, we invoke a gauge transformation (first 
proposed in Ref. \2% in the Usadel equation H3.30|l : 



so that 



where the matrix field 

'K+ K. 



K = 



(4.1) 



(4.2) 



K 



will be chosen such that almost all of the contributions 
of <f> from frequencies u> > 1 /r, are eliminated4£*21 The 
corrections to g can then be found by perturbation the- 
ory. 

We substitute Eqs. (|3~32'|) and gTj into Eq. (pOnjl and 

obtain 



(4.3) 



.g;e tK f n e 



1 + — [\9> e Ue }-2 



We will look for a perturbative solution to Eq. I|4.3I) in 
the first loop approximation; to do so it suffices to retain 
terms that are at most quadratic in the field K on the 
right hand side. At zeroth order in </>, the Green function 
has the form 



fJ 







t 2 ) g K (ti. 
-6(h - 



(4.4) 



K 



We will require that this form is preserved even in the 
first and second orders, i.e. the corrections to the spec- 
trum (given by g R and g A ) are indeed eliminated by the 
gauge transformation. Considering the R, A and Z com- 
ponents of Eq. (|4.3|) in linear order in K , we see that 
this condition is satisfied if 



| z) </)=<<_(/). 



(4.5) 



K- = -C 9 <f>- 



(4.6) 



Here the bar indicates complex conjugation, and the 0- 
dimensional diffuson C a is given by 



£»H = — 1 



(4.7) 



-iu> + 1/r ' 

with 1/t the escape rate Eq. (I2.1U|) in units where 6± = tt 



1 \ - T n _ g L + gR 

-r ^ 2 ~ 2 ' 



and we introduced the operator notation 



( tl )= dta&fa-taW-fa). 



(4.8) 



(4.9) 



The diffuson C 9 {t) gives the classical probability for an 
electron introduced to the dot at t = to remain in it 
at time t. Note that in the present case the diffuson is 
much simpler than in higher dimensional systems, where 
it depends also on position in space and direction of mo- 
mentum on the Fermi surface. The same is true for the 
fields <f> and K. 

To solve the Kcldysh component of the gauge trans- 
formed Usadel equation l|4.3|l . we will look for K + in the 
form 



[C 9 



2K_ 



K- = {idt)~ l MK_, (4.10) 



where the operator M(t±,t2), will be chosen to simplify 
further expansion of g K . Substituting Eq. I|4.1Q(I into the 
Keldysh component of Eq. 14. 3D . and expanding the right 
hand side to second order in K we obtain: 



dg* 



Of 



+iJ2^[K + ;f«-g K }+iQK_-2i 

n 

-Et{[^[^;/"]] 

n 

+ (1 - T n ) [K+tf K^g K - g K K^f«K+] 
+T n [K + f «K_f « - f«K_f«K + ] }. 

Here we do not display second order terms that vanish af- 
ter averaging over the field if>, e.g. K-K- (see Eq. (14. 18(1 ) 
and K + (t 3 )K-(t 3 ). The latter product vanishes due to 
the analytic properties of the retarded propagator (see 
Eqs. (|4~T%1 and J03JI). Defining 6g K to be the linear 



(4.11) 



K-;g 



K 
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order correction to g , we find: 

l + l)sg« = ij:^[K +] f«- 9 «}, (4.12a) 



2? 



(4.12b) 



where the if's are understood as operators: K(ti,t 2 ) = 
K(ti)5(t\ — t 2 ). The operator Q(ti,t 2 ;t 3 ) acts on as 



[QA-_](ti,t 2 )= j dt 3 Q{tx,t 2 ;t 3 )K^{t 3 ), (4-13) 

and is related to products of the Keldysh Green function 
of the dot and the leads: 

Q(h,t 2 ;t 3 ) = (4.14) 
\Y.{ T -[fn{tiM)f«{hM) + g K (tiM)g K {hM)] 

n 

+T n (l-T n ) [(g K - f«)(h,t 3 )(f« ~ g K )(t 3 ,t 2 )]} . 

The Keldysh Green functions are singular at coinciding 
times: 



g K {t l ,t 2 )\t 2 -,t 1 = - 



2i 



7T(tl - t 2 ) 



+ regular, (4-15) 



(similarly for ) and therefore the products in the sec- 
ond line of Eq. 1|4.14(1 are understood as 

g K \hM)g K '(h,h) = 

\ g K (t 1 ,t 3 + <7iO)g K (t 3 + <7iO ) t 2 ), ( 4 ' 16 ) 

1 <r=±l 

and similarly for • . We will see later that the 
operator Q has the meaning of the fluctuations of the 
charge going into and out of the dot. The third line in 
Eq. (|4.14|) vanishes in equilibrium and represents the non- 
equilibrium shot-noise, whereas the second line is present 
even in equilibrium and has the meaning of Nyquist noise. 

As the right hand side of Eq. H4.12al) vanishes for 
t\ = t 2 , we see that Sg+(ti,ti) = 0. This means that 
5g+ does not contribute to physical quantities of the dot 
itself (e.g. total charge or total energy in the dot). Fur- 
thermore, taking the limit t\ — > t 2 in the right hand side 
of Eq. Q4.12bjl and using Eq. (|4.15() . it is seen that the 
same condition can be imposed on 5g^_ if we choose the 
operator M to be 



M{h,t 2 ) 



--Q(h,h;t 2 ). 



(4.17) 



We note how the structure of the Q and M operators is 
similar to that of Eq. (5.30) of Ref. [H for the case of 
disordered metals, the difference being that in that ref- 
erence the Green functions are characterized by different 
directions instead of different channels. 



The next step is to average Eq. (|4.11(l over the 
Hubbard-Stratonovitch fields. For this we need the prop- 
agators for the fields K ± defined as 



(K + (t 1 )K + (t 2 )) !l> =-K K (t 1 ,t 2 ), 
(A' + (ti)A-_(t 2 )) =iA- R (ti,t 2 ) J 



(AL(t 1 )A' + (t 2 )) =-l^(t 1 ,t 2 ), 
(X'_(ti)A'_(t2))^=0. 



(4.18) 



Their relations to the propagators H3.1fi(l and their ex- 
plicit expressions will be given in subsection IIV A2l 

After averaging over the fluctuating fields and drop- 
ping from now on the superscript K in the averaged 
Green functions (g = (g K }<p), Eq. (|4.11|) takes the form 



dg 

at 



E IT (/» " 9) + St! + St a , 



(4.19) 



where the collision integrals Sti and St 2 arise from aver- 
aging terms of first and second order in K± respectively. 
To average the second order terms (given by the last three 
lines of Eq. 1)4. llf) ). we replace g K with its average g - 
keeping the deviation from average will result in higher 
order terms beyond our approximation- and average the 
K± pairs using the propagators (|4. . obtaining 



St 2 = St 



St 

i 



(4.20a) 
(4.20b) 



o J2[ T nK K {h,t 2 )f n (h,t 2 ) 

n 

l -T n (1 - T n ) J dt 3 (K* 3 - K£ 2 ) (A" 3 g 3>2 + g h3 % 2 ) 



Kl 2 )fl 3 f- 



n fn, 
3,2 



-T^Jdh (Jf« 3 - 

&*{tx,ti) = ~^T n (\-T n ) (4.20c) 

n 

x J dt 3 (A" 3 - «/;., //•,,) (a , /; :1 + Kl 2 ) , 

where 

K K (h,t 2 ) = 2K K (hM) - K K {t u t x ) - K K (t 2 ,t 2 ) . 

(4-21) 

As for the first order terms (given by the second line of 
Eq. 14. llll ). we write g K = g + Sg K , where 5g K given 
by Eqs. (|4.12a(l - l|4.12b|l is the correction linear in K±. 
Keeping terms of second order in K±, we obtain 



st 1= (st_) + (st+; 

St_=i [SQ] K_ - 2i 
St+=-- [K + ;Sg K ] 



K 



(4.22a) 
(4.22b) 

(4.22c) 
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where SQ is the first variation of Q [Eq. Q4.14J1 ] with 
respect to g 



K. 



SQ 
Sg 



K 



Sg 



K 



(4.22(1) 



Using the solutions of Eqs. H4.12ajH4.12bf! f° r Sg K , one 
can explicitly average St± using Eq. (|4.18|l . but as it will 
be shown in the next subsection the collision integral Sti 
does not contribute to physical quantities and we will not 
need its explicit form. 

This formally concludes the derivation of the equation 
for the Keldysh Green function in the time domain. How- 
ever we need to calculate explicitly the propagators intro- 
duced in Eq. I|4.18|l : this is done in the next subsection. 

2. Propagators and Collective Excitations 

The propagators defined in Eq. I|4.18(l can be expressed 
in terms of the propagators V, Eq. (|3.16|1 . thanks to the 
relations between the fields K± and (f>±; see Eqs. (14 .6[) 
and H4.10fl . We obtain 

r r = _ C g T) R C g j (4.23) 
K K = C 9 V%C 9 

+2i [L g dl x MK A - K R Md 1 T l C 9 ] . (4.24) 

Therefore we need to evaluate the T> propagators, given 
by Eq. (pTTTjl and Eq. (fTMjl. W e note that the Green 
functions appearing in Eq. (|3.34|l are the ones before the 
gauge transformation Eq. (|4.1() . so that to linear order 
in (j), the Sg's obtained from Eqs. I|4.12a|l - I|4.12b(l have to 
be modified according to 

8g K -> Sg K - i[K+-g K ] - 2K.5(t 1 - t 2 ), 

5g z -► 2KS(ti - t 2 ). (4.25) 

We also note that, thanks to Eq. (|4.15|) : 



lim -i[K + ;g ] = d t K^ 



(4.26) 



Using this result together with Eqs. (|4. — (|4. 1U|) and the 
property 

Sg K (t 1 ,t 1 ) = 0, (4.27) 
from Eq. i|3.34[> we find 

n£(ii, t 2 ) = ~[&{ti - ta) - d tl £ 9 (tx, ta)] , (4.28) 
or after the Fourier transformation (|4.7|l 

r 7T — IUJ + 1/T 

The result for the Keldysh component is 

v 7T 



(4.29) 



(4.30) 



Substituting Eq. (|4.29() into the retarded component of 
Eq. <|3.17j) we obtain 



2>£(w) = 2£ c 



— iuj + 



(4.31) 



where we restored dimensionful units (the difference be- 
ing the mean level spacing appearing in the denominator 
instead of 7r). Then we can calculate the K R propagator 
in Eq. I|4.23|l : it can be written in terms of the diffuson 
("ghosts") propagator C 9 and the following propagator 
C p for the collective excitations in the charge channel 



as 



Si 1 

2 — iuj 



C 9 {lo) 



1 



~6\ 



(4.32) 



(4.33) 



To find K K we use the Keldysh part of Eq. ll3~T7|) : 

Vf=VflL«V$, (4.34) 
together with Eqs. (|4.24|) and (|4.3U|I and obtain 



K 



K 



l 2d t 

- [ 1 + 



4£, 



£W 9 + AT S £ S 
C P N P +N p £, p 



(4.35) 



1 



written in the time domain and in the operator notation. 
Here the bosonic "density matrices" N a , a — g,p, are 
defined through 

(C 9 )- 1 ^ 9 +Af 9 {£ 9 )- 1 = 2M, 

{C^N" +Af p (£ p )- 1 = 2(1 + F p )M, (4.36) 

i.e. they are required to satisfy the kinetic equations 



-2(1 + F a ) [ -N a 



M 



with 



F 9 = , F" 



AE C 
Si 



(4.37) 



(4.38) 



While in general we will not solve these kinetic equa- 
tions, their exsistance is needed to obtain the conserva- 
tion law for the energy. To convince the reader of the 
bosonic nature of the collective excitations, let us briefly 
consider the thermodynamic equilibrium 17(e) = / n (e) = 
2tanh(e/2T). In this case, after Fourier transforming 
Eq. (l4*T7|l we find 



1 , w 
-wcoth— , 

T Zl 
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and introducing the distribution functions N a (t,to) by 

Af a (h,t 2 ) = J ^e- %uj{tl - t2) ,Lo\2N a (Np.,u) + 1 

(4.39) 

we arrive at the followig solution of Eq. (|4.37|) : 



(4.40) 



i.e. the Planck distribution. 



3. Conservation Laws and Currents 

The validity of the conservation laws in the kinetic 
equation approach is related to certain properties of the 
collision integral. For example in a closed system the 
charge conservation law follows from the vanishing of the 
collision integral in the limit t% — > ti 16 In the present 
case however we are dealing with an open 0-dimensional 
system and therefore the conservation law for any phys- 
ical quantity like the charge q(t) in the dot should have 
the form: 



dq 



0, 



(4.41) 



where I = Il+Ir represents the total charge flux leaving 
the dot. 

To obtain the charge conservation law, we rely on 
the following general properties of the collision integrals 
(I4~2tl - (l4~m 



lim St =0, 

*2->tl 

lim St+ = 0. 

t 2 ->ti 



(4.42) 



The last property is a direct consequence of the definition 
(|4.22c[) together with Eq. The proof of the first 

property can be obtained by following the steps described 
in Appendix E of Ref. 0. We also notice that for the 
stationary state [see Appendix |A"] : 



lim St_ = 0. 

ta->ti 



The charge in the dot is given by 



30 



Z *2— >tl 



(4.43) 



(4.44) 



and taking the same limit of both sides of the kinetic 
equation 1)4. 19[) we obtain Eq. (|4.41|l with: 



I = Il + Ir + 1- 

de [T„, 
2^ 



(f n (e)-g(e))+I^(e) 



I- = e lim St_, 



(4.45a) 
, (4.45b) 

(4.45c) 



where a = L, R and 



x [/„(e - u)g(e) - g(e - w)/„(e)] 



(4.45d) 



The contribution /_ vanishes in the steady state, thanks 
to Eq. (|4.43(l . and will not be given any further con- 
sideration. We mention that the above expressions for 
the electric current can also be obtained by applying the 
gauge transformation Eq. 1)4. lfl to the current formula 

Eq. arm . 

We now turn to the energy conservation law, whose 
validity is based on the properties 



hm - (d tl - d t2 )St m = dtMti) +E ^(ti), (4.46a) 

n— 1 

lim (^ 1 -9 t2 )St+ = 0, (4.46b) 

h '2— *tl 

where 



-t dt K K d t -\ Jdt 3 [d t K R ] 13 g h3 g 3A 

~^T n (l - T n ) J dt 3 [d t K R ] 1 3 [«?(/„ - g) + (/„ - g)g] , 

(4.47) 

or after the Wigner transform 
duj 



I 



4tt / 2tt K ' 



(4.48) 



T n (l - T n ) f dedu 
4^2 



■uImK R (ij)g(e-u)(f n (e)-g(e)). 



As before, the last property can be proved straightfor- 
wardly, while the derivation of the first one is delineated 
in Appendix E of Ref. 1 161 In order to write it in the above 
form we used the kinetic equations (|4.37(l and defined Ub 

as 



Ub = Up - Ug, 

u a = ^C a Af a , a = p,g. 



(4.49) 



Below we will identify it as the contribution of collec- 
tive excitations to the dot's energy. The separation in 
Eq. I|4.46a|) of a energy contribution from a energy cur- 
rent one may seem arbitrary; however there are two inde- 
pendent tests of the validity of Eqs. (|4.48(l - (|4.49(l . First, 
since the leads' electrons are non-interacting, we can pro- 
ceed similarly to the derivation of Eq. I|3.31|l for the elec- 
tric current and find that the energy current is obtained 
by replacing e — > i(d tl — d t2 )/2 in that equation; sub- 
stituting in the resulting formula Eqs. I|3.32|l and Q4.1[l . 
expanding to second order in K± and averaging over the 
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fluctuating field, we again arrive at the results presented 
below, Eqs. I|4.54|l . Second, the definition l|4.49|l of the 
collective excitations' energy is in agreement with the re- 
sult of the thermodynamic calculation for such quantity, 
see Appendix [5] Finally, we note that the property sim- 
ilar to Eq. (|4.43(l holds in the steady state: 



lim (d tl - d t2 ) St_ = 0. 

*2— >tl 



(4.50) 

The electrons' contribution to the dot's energy u is 



u e (t 1 ) = -- lim (d tl -d t3 )Tr ll g{t 1 ,h), (4.51) 

4 12— >tl 



and from Eq. 14.19f) we arrive at 

dtutot + I 6 = 0, 



(4.52) 



where the total energy utot is the sum of the contributions 
of the electrons and collective excitations: 



Utot — u e + u b 
The energy current I £ is: 
I 



Il+I £ R + I- 



de 



2ir 



T, 



-f(fn(e)-g(e)) + I* l (e) 



It = - lim (d tl -d t2 )St. 



(4.53) 
(4.54a) 

je.in 
" l n ' 

(4.54b) 
(4.54c) 



where a = L,R, 1^ is defined in Eq. (|4.45dl) and I% m in 
Eq. (|4.48|l . Again, since It vanishes in the steady state 
- see Eq. (|4.50|) - we will not discuss it anymore. 



4- Final Form of the Kinetic Equation 



and 



duj 1 

2-7T LO 



(4.57) 



£"Mf«( e ,c) + f«( £ , w )£«M 



T n (l-T n ) 



Re£ a (uj) n n (e - oj)(n n (e) - n(e)) 



where the propagators C a , a = g,p, are defined in 
Eqs. I|4.7|) and l|4.32[l . and we introduced the combina- 
tion of distribution functions: 



T£(e, lo) =(N a (u) + l)n n (e)(l - n„(e - w)) 
-N a (io)(l-n n (e))h n (e-uj). 



(4.58) 



The bosonic distribution functions N a (u), a — g, p, were 
defined in Eq. I|4.39|l and they satisfy the kinetic equation 



8N a (t 



at 



1 + F° 



'T«(e )W )+T«( e>W )" 



+ (1 - T„) [n(e) - u„(e)] [n(e - u>) - n„(e - w)] L 



(4.59) 

which follows from Eq. (|4.37|l . Here, T(e, u>) is given by 
Eq. I|4.58J) after the replacement h n — > n. The combi- 
nation 14.58|l can be obtained by the standard argument 
for the creation and annihilation of the one-particle ex- 
citations in the dot-lead system. On the other hand the 
last line in Eq. (|4.57(l may be understood as the renor- 
malization of the scattering coefficients T n of the non- 
interacting collision integral Eq. I|4.56|l . due to interaction 
with the self-consistent potential in the dot. 



Closing this section, we present for completeness the 
final form of the kinetic equation for the dot's distribu- 
tion function n(e) — \ — \g{e). The collision integral 
is a functional of n(e), the leads' distribution functions 
«n(e) = \ ~ i/n( e ) an d the bosonic distribution functions 
N a (ui) defined in Eq. g^: 

^± = §t{n,h n ,N a }, (4.55) 
ot 

St = St T {n, h n } + §t P {n, n n , N p } - St 9 {n, n n , N 9 }. 

In the second line, we separated in St three physically 
distinct contributions. The first term on the right hand 
side describes the relaxation of n due to tunneling into 
and out of the contacts - this mechanism is present even 
for non- interacting electrons. The two other terms de- 
scribe the interaction effects. Explicitly, they are given 
by 

St T (e) = -£^[n(e)-n n (e)], (4.56) 

n 



B. Triplet Channel 

The case of interaction in the triplet channel can 
be treated simlarly to the singlet channel, the main 
difference being that the Hubbard-Stratonovich field, 
Eq. i|3.37fl . is now a vector. Therefore the phase fac- 
tors in the gauge transformation must also possess this 
structure and in Eq. I|4.1|) we substitute 

K(t f ) K(ti)®a. (4.60) 

This transformation does not commute with the Zeeman 
energy term in Eq. (|3.42l) ; after the gauge transformation 

and expanding up to second order in K this term becomes 

l ~E* z [a z :g] - ±E% ( [a z ;g] + iK a [[a a ; a z ] ;g}), (4.61) 

where we sum over the repeated index a = x,y, z. Note 
that the second order terms vanish identically. 
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The linear order equations for K± can be decoupled by 
using the following basis for matrices in the spin space: 



A = ^Tr s (A)l s +A m a r 



A m = ^Tr s (a- m A) 



c = o z , o =-={a x ± ia 9 ) . 
V 2 



(4.62) 



These matrices obey the following commutation rela- 
tions: 



and the density matrices a = g, s, are defined simi- 
larly to Eqs. P5) . 

We could also repeat (with necessary modifications) 
the calculations of Sec. IIV Al and find the final form of 
the fermionic collision integral. The kinetic equation 
Eq. I|4.55J) . is modified by adding the following collision 
integrals to the right hand side: 



E 

m=0,±l 



(st m st m ^j 



(4.72) 



[a ; a* 1 ] = ±2cj ±x , [<t +1 ; a' 1 } = 2a° , (4.63) 

and all other commutators vanishi^i Then the solution 
for K m , m = 0,±1 is 



K" 



-CLh m 



(4.64) 



(no summation over to) , where C 9 m is obtained by shifting 
the frequency of C 9 in Eq. I|4.7|) : 



(4.65) 



We can similarly obtain the expressions for K™ . In 
the new basis, the formula for the polarization opera- 
tors n mj _„ is given by Eq. (|3.39[) with (i, j) replaced by 
(to, n), and we similarly obtain 



fri? 



and using Eq. 



imE% 



J, 



L "J mn v ' 9 



7T — i(uj — mE* z ) + - 
-i (u> -E* z ) + ± 



1 u m 7 —n i 



(4.66) 



2 + (I+ztoS*) 

(4.67) 

for the interaction propagators. Using these expressions 
in the zero frequency limit, together with Eqs. I|3.41|) - 
(|3.43|> we obtain the expression Eq. I|2.11[) for E* z . Finally 
for the retarded propagator we obtain K R n — K R 6 m .- n , 
with 



where 



Si 1 

2 — iui 



C° m (u)-(1 + F')£' m {w) 



-i{Lu-mE z ) + (l + Fs)±' 



and 



Si 



Similarly, K* n = K*5 m ,- n where 



2 d t 
(1 + F S ) 



£9 Kf9 i ATS £9 
J -'m JX 'm ~ - ,v m'—m 



(4.68) 



(4.69) 



(4.70) 



(4.71) 



£ s W s 4- J\f s C s 



1 

ft' 



where St m are given by Eq. (|4.57(l after replacing £ Q with 
and iV Q with N*. The AT£, a = g,s, are distribu- 
tion functions for bosons with unit spin and are defined in 
terms of the density matrices Af m as in Eq. (|4.39|) . They 
satisfy the kinetic equation Eq. I|4.59() after the replace- 
ment N a — * N m . It also follows that the expressions for 
the currents are obtained by replacing 



K 1 



in Eqs. (|4"43|) and l(33i|) . 



m=0.±l 



(4.73) 



C. Cooper Channel 

For the interaction in the Cooper channel, we resort to 
a perturbative approach which, in contrast to the charge 
and triplet channels, does not start with a gauge trans- 
formation. Instead, we use the constraint (|3.27|) . to write 
the retarded and advanced Green functions as: 



1 - \T R -T l 
T R 



T R 



A + \t R -T R /N 



\T A -T A 
T A 



1 



I 'pA -rA 
2- r ' / N 



(4.74) 



(4.75) 



F K = +~{F R -9 K + 9 K -F A ), 



(4.76) 



and the Keldysh anomalous Green functions as: 
1 

2 ( 
1 

'2 1 

where we use the notation Eq. I|3.46|l . In the normal 
phase the anomalous Green functions have only fluctu- 
ating parts proportional to A, and hence the above ex- 
pressions are valid up to second order in the fluctuating 
field. Using these formulas together with Eq. I|3.32|l . and 
defining 

N ch 

Sti = -i[A;g], St 2 = (4.77) 

n=l 

we can write the upper left GN component of the Keldysh 
part of the collision integral St 2 as 



(I, 



fJ 



St^ ~\~ St r 



- MT 



(4.78) 
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where we have identified the Maki-Thompson and elastic 
parts as 

- MT 1 



St, 



St. 



l -T n {T n -l)[T R T R {f: 



■(/, 



K 



i K \T A T A 



(4.79) 



(4.80) 



We note that there is no term analogous to the Asla- 
mazov-Larkin 28 contribution this is due to the indepen- 
dence of the field A from the spatial coordinate. Simi- 
larly, the substitution of Eqs. I4.74[) - I4.76[> into Sti gives: 



(stoft = -%a+t 



K 



iF K A* 



iA_T A - iJ 7 A* 



l -A + T R g K ~ l -A + g K T A 



F R g K A% 



--g K T A A% 



2 

iA_f A 



iT R A*_ 



(4.81) 



In order to average over the fluctuating field A, we 
need to solve the equations of motion for the JF's. Using 
Eq. H3.53fl for the T R component, we obtain 



9t, - d,. 



F R + l -E z {a*;F R } 



tJ 2 > (4.82) 

+2iA+S(t 1 - t 2 )+ig K A_, 



where 



1 

T* 



gh 
2 ' 



(4.83) 



and gh was defined in Eq. (|1.2fl . This modification of the 
escape rate is due to the last term on the left hand side 
of Eq. (|3.53|) , which describes the breaking of time rever- 
sal invariance for the anomalous Green functions by the 
orbital magnetic field. It is clear from this equation that 
the .F's are diagonal in spin space, and if we decompose 
them according to 



(4.84) 



the solution to Eq. (|4.82|) can be expressed in terms of 
the cooperon C as: 



T R (h 

+i J dt 3 C± 
where 



t 2 ) = *C±i 

h-t 2 



C±(t) = 



2 

*a) A_ 
de 



h + t 2 



(4.85) 



t 3 )g K (2t 3 ), 



(4.86) 



The T R is obtained similarly and is given by Eq. (|4.85() 
after replacing A± by tA±, and g K by g K . For the 
advanced components we can use Eq. (|3.47|) to write: 



F£(ti-t 2 )=F R {t 2 -t 1 ); 



(4.87) 



the same relation holds for the JF's. From Eq. 13. 4711 we 
also have 



-K 



(t 1 ,t 2 )=g K {t 2 ,t 1 ), 



(4.88) 



and similarly for /„ and /„ (here and for the rest of this 
section we drop the superscript K from the dot and leads 
Keldysh Green functions). 

Using these solutions for .F's and the propagators 
(|3.54|) to average over the fluctuating fields, we arrive 
at the following expressions for the singlet part of the 
Fourier transformed collision integral St2: 

-Tr^St^e)) = \ T nJ ^(/»(« " e) - - e)) 

x Y, \C m \ 2 (2e~u } )J(e,o J )\V R (u;)\ 2 , (4.89) 
m=±l 



iTr s (sC'(e)) 



±T n (l-T n )[^(f n (e)-g(e)) 



^ m=±l 



u)J{u-e,u>)\V R {u)\ 



+g(co-e) J2 I™[C? n (Lo-2e)V R (Lu)] 

m=±l 

with the kernel J(e, uj) given by 

J(e, u) = »Il£(w) + g{e)ImIL%(uj) . 



(4.90) 



(4.91) 



To write the collision integral in this form we used the 
identity 



\mK(u) + g(e)lmIl R (u)]\D R (u)\ 2 
= W*[(u)+g(e)lmVZ(u), 



(4.92) 



which follows from the Dyson equation H3.55J1 for the 
propagators. At uj much larger than T or applied voltage 
eV, iT>^(uj) = — 2sgn(w)ImZ>§(u;), and the contribution 
from such large w's in Eq. H4.89JI vanishes. This is a direct 
manifestation of the inelastic nature of Maki-Thompson 
processes. 

By comparing Eqs. (|3.31|) and l|4.78|l with the defini- 
tion of St2 [Eq. H4.77fl ] we find that the current in each 
channel is given by 



e f de 

In = 2 Trs / Y 



T n el 

jiz-i-jj + st,, 



St,, 



■ MT 



(4.93) 



with the expressions (|4.89|) - (|4.90l) for the collision inte- 
grals. We will use this result in Sec. lVGI to calculate the 
corresponding interaction correction to the conductance. 

By repeating the above steps for the collision integral 
Sti we find: 



1 



Tr s (St!(e)) 



1 E 

2 ^ 

m=±l 



du) 



Re(C m (2e-w))£(e,o;) 

(4.94) 
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where the kernel /C(e,w) is 

K(e,u>) =iV%(u;)[g(e) + g(co - e)] 

+ lmVi(cu)[g(e)g(u J -e)+A] 



(4.95) 



Having derived the collision integrals, we now find the 
explicit form of the propagators Eq. (|3.54|) . and polariza- 
tion operators Eq. (|3.56|) . Using Eqs. (I4.76[) and (|4.85|) 
we obtain 

Ha(w) = y /I 1 E C ro (2e) 5 (|- e ), (4.96a) 

•* 77 m=±l 

nlH = 5 E C ro (2e) 5 (| +e), (4.96b) 



de 
2^ 



E C m (2e)[ 3 (|-e).g(^ 



rn — ±1 



+ e +4 



(4.97) 

In Eq. 14.97JI . the last term in the brackets does not follow 
from the Usadel equation. It arises to compensate for the 
incorrect order of summation over orbital states m and 
integration over e in Eq. (|3.56|) . similarly to the anomaly 
in charge channel [cf. the discussion after Eq. QH.34jl ]. 
This ultraviolet factor can be obtained e.g. by explicit 
calculation in the original model. A more compact way, 
however, is to restore the correct term by requiring the 
fluctuation dissipation theorem: 



ntH = coth(^) [ngH-iElM] 



(4.98) 



to hold for the equilibrium "distribution function" g(e) = 
2tanh(e/2T). Furthermore the logarithmic divergence in 
Eq. (|4.96|l is cut off by |e| < Et as it determines the 
validity of RMT^ 

In equilibrium, with g(e) given by the right hand side of 
Eq. 1)3. 8JI at temperature T and [i = 0, we can calculate 
explicitly by closing the integration contour in the 
upper half plane, obtaining a sum over the residues at 
the poles of tanh(e/2T). The result is 



n e »= E <(«) = -l E Err r 



m=±l 



27r ^— ' ^— ' n + 7j + z Tl 

m— ±1 n— ^ 



where 



-i{uj - mE* z ) + ± 
4^T 



(4.99) 



(4.100) 



The upper cutoff N is approximately given by Et/2ttT . 
We note that the coupling constant J c is also defined at 
the scale of Et ■ We can express this sum in terms of the 
digamma function i/j(°\x) 

- 2.11* H = ^ (\ + ^+ J - ^ (\ + *, 



In 



\2ttT 



V^ (0) ( ~ + 2, 



(4.101) 



Defining the "critical temperature" 

2-yEr 



we can rewrite IT5; as 



(4.102) 



2ttII* ( W ) = n J- 1 + t + ^ Q + z m ) - V (0) (0 , 

(4.103) 

where t is the reduced temperature 



, T 
t = In— . 

T c 

Therefore, in equilibrium we find 
KM = -TTT- 



(4.104) 



<+li: m=± i^ o) (i+^)-^ o) (i)' 

(4.105) 

Furthermore, in equilibrium we have the fluctuation- 
dissipation relation: 

KM) = coth (^) \D%{w) - 2? e »] , (4.106) 

as can be verified using Eqs. I|4.96[) . I|4.97[) and l|3.55[) . 

We can define the bosonic distribution function of the 
fluctuating cooper pairs by generalizing this relation to 
the non-equlibrium case: 

Vliuj) = (2iV» + 1) - ZtfM] , (4.107) 

such that in equilibrium N c (lo) is given by the Planck 
distribution Eq. (|4.40|l . 

The polarization operators must be recalculated in the 
non-equilibrium case, as they depend on the distribution 
functions of the electrons. We do this for the DC case in 
Sec. |VT?] 

In closing this section we give the final form of the 
kinetic equation for interaction in the Cooper channel. 
Introducing n(e) = \ - \g(e) and n„(e) = \ - \f n (e) 
for distribution functions of the electrons in the dot and 
the leads respectively, and using Eqs. (|4.89|) - (|4.94|) and 
(|4.107|l we arrive at 



dn{e) 
dt 



T 

= - E ~Y ["( e ) ~~ "«( e )] + St e({"-, n n , N c } 

n 

+St M T{n, n„, N c } + Sti{n, h n , N c }. 

(4.108) 

The collision integrals are given by 
St eI = -5^T n (l-r n )[n(e)-n n (e)] / '^M(e,w), 

n J n 

(4.109) 

where 

X(e,c) = [iV c H+n(e)] E Ref* (2e - w)Im2>2(a;) 

m=±l 

-\n{u> - e) E Im ~ , 



rn — ± 1 



(4.110) 
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and 



n m=±l 

x [* n ( e> w)-*(e,w)] , (4.111) 



— ^ ReC ro (2c-w)Im2?f(a;)*(e ) 

77 m=±l 



Here we defined 

*(e,w) = - (N C H + l)n(e)ra(a; - e) 

+ 7V c H(l-n(e))(l-n(c- e )), 



(4.112) 



(4.113) 



and ty n is obtained by replacing n[u) — e) with n n (u> — e) 
in Eq. If4.113fl . The bosonic distribution function N c (uj) 
is in turn given in terms of the fermionic distribution 
function n(e) through 



2N c (lu) + 1 = 



(4.114) 



and the expressions H4.96fl ~ H4.97fl f° r the polarization op- 
erators. 



V. CALCULATIONS OF THE TRANSPORT 
COEFFICIENTS 

In this section we explicitly calculate the interaction 
corrections to the electrical and thermal conductances 
by solving the kinetic equations. We first consider the 
singlet and triplet channel corrections to the electrical 
conductance and then to the thermal conductance. Fi- 
nally, we calculate the Cooper channel correction to the 
electrical conductance. 



A. Electrical conductance: singlet and triplet 
channel 



In the presence of a bias voltage, the energy levels in 
the leads are shifted by the voltage times the electron's 
charge, so that the nonequilibrium "distribution func- 
tions" in the leads, Eq. 1)3. 8fl . become 



( e ) = 2 tanh 



e-eV, 



R(L) 



2T 



(5.1) 



with Vr(£) = ±V/2. For the future use we introduce 

A/(e) = f L (e) - Me). (5.2) 

In the steady state we have II = —Ir = I, so we can 
write 



I = xI L + (x - 



R, 



(5.3) 



for any constant x. In terms of the current in each chan- 
nel 



(5.4) 



n=l 



where we choose the constants to be 



gR 



A,, 



Sl+Sr 



1 < n < N L , 



N L + 1< n < N ch , 

gL+gR ^ — — Cn ' 



(5.5) 



with gi, gR defined in Eq. I|2.3fl . This choice will simplify 
calculations as we make use of the identity J~) A n T n = 0. 

Let us calculate the singlet channel correction first. 
Using Eq. (f4.45b() we have 



A',, 



de 



(5.6) 



with J®' given in Eq. 14.45dfl . By shifting the integration 
variable e by ui in the first term in square brackets in 
Eq. 14.45dl) . we can rewrite I as: 



I=e 



de_ 
2^ 



T, 



£ A „ ^(/„(e)- 5 (e)) 



1, 



+ -T n (l-T n )f n (e)K(cj)g(e-u) 



(5.7) 



where 

K(w) = -ImK R (u) = Im[(£ 9 ) 2 {uj)V r (uj)} , (5.8) 

and we used Eq. (|4.23l) in the last identity; in this form 
we can recognize that our expression for the current has 
the structure similar to the one found in Ref. 8 for higher 
dimensional systems. However, for practical purposes, we 
will use the following identity: 



K(ui) = —Im[£ 9 (u;)-£ p (u 



(5.9) 



where C g and C p were defined in Eqs. I|4.7fl and (|4.32fi 
respectively. 

We note that since J2 n ^nT n g = 0, the interaction cor- 
rection to g does not enter the expression for the current 
to first order in the interaction propagator. So for the 
purpose of calculating the DC current, we can calculate 
g to zeroth order by equating to zero the first term on 
the right hand side of Eq. I|4.19fl and we find: 



9(e) = 



gL/L(e) + g fl /fl(e) 

gL + gR 



(5.10) 



with gL,g_R defined in Eq. I|2.3fl . 

We write the differential conductance as 



G 



dl 
~d~V 



Go + AG, 



(5.11) 
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where Go and AG are the classical conductance and the 
interaction correction, originating from the first and sec- 
ond term in curly brackets in Eq. (|5.7[) respectively. Sub- 
stituting Eq. H5.10fl into Eq. 1)5. 7f) and using the definition 
(15.111) we obtain 



Go 



e 2 gLgR d I" de_ 
2 gL + gRdVJ 2tt 



A/(e) 



iLgR 



7T gL + gli 



(5.12) 



and 



AG = 



g|h_L + g|h fl e 2 

(gL 



IR) 2 4 



dVJ 2ttJ 2tt w v m ; ' 



(5.13) 



where the form factors h^ and h# are defined in Eq. I|2.7|l . 
and 

V = eV. 

We define the shorthand notation 



aKb= j^J ^a{e)K(u;)b{e-u), 



(5.14) 



and notice that since K is an odd function of u), we can 
further simplify Eq. I|5.13[l into 



, r e 2 h L g 2 R + h R g 2 L 7T 

7T (gL+gfl) 4 



(5.15) 



where here and below a prime denotes derivation with 
respect to the energy variable. The integration over e 
can be readily performed and after a partial integration 
over to we arrive at 



2T 



2tt v ; 2T 



coth ■ 



r 



2T 



(5.16) 



Using the identity (|5.9|) together with the definitions 
Eqs. (|4.7|) and (|4.32|) . we integrate over cu by closing the 
contour in the upper half plane, obtaining a sum over the 
residues of the poles of coth (w/2T), which is expressible 
in terms of polygamma functions. The result is: 



Re 





-( 




y 2htT J 





1-ktT 



(5.17) 



where 



(5.18) 



and i(}^(z) is the i— th derivative of the digamma func- 
tion. Substitution into Eq. (|5.15|l yields Eq. I|2.5f) for the 
singlet channel correction to the electrical conductance 
(where we have also restored dimensionful units). 



As for the triplet channel, the result is found by using 
the replacement in Eq. 1|4.73[) together with the defini- 
tions (SEjg) and g203>- Then Eq. becomes: 



/t + im_E; 



-Im[£^(o;)-£^H] 



m=0,±l ""' ' ' """ z 

(5.19) 

and repeating the above steps we obtain Eq. H2.tjfl for the 
triplet channel correction AG S . 



B. Thermal conductance 

The thermal conductance is found by assuming that 
there is a temperature difference ST between the right 
and left leads, so that the distribution functions, cf. 
Eq. Q3.8p. are given by 



Ir(l) (e) = 2 tanh 



2T±5T 



(5.20) 



where + (— ) signs corresponds to right (left) lead. At 
lowest order, g is still given by the relation l|5.10[) . 
We write the linear thermal conductance as 



k = t^; = k + Ak c1 + Ak;„, 
61 



(5.21) 



where ko is the classical thermal conductance and An e i 
and AKi n are the contributions of the elastic and inelastic 
processes described respectivley by the collision integrals 
in Eqs. H4.2()c|) and l)4.20bjl . The calculation is similar to 
the one in the previous section: we can write the steady 
state energy current as A n I^ and from Eqs. 
and H4.54f> we obtain 



I E — ^0 + lei + lin 



(5.22a) 

T o = o E A « T » / Y £(/n(e) " 5(£)) ' (5 - 22b) 

N h 



n=l 



xe[/ n (e - w)ff(e) - <?(e - w)/»(e)] , (5.22c) 
1 . , , f de f duj d, . 

^=-i6Ew-T n) y-y-^v) 

n=l 17 17 

x [f(e - w) (/n - 9) (e) + (/n - ff) (e - w)s(e)] • 

(5.22d) 

We note that the part of If proportional to T n 's [see 
Eq. (|4.48(l j. vanishes in the DC case after summing over 
the channels due to the identity ^ n A„T„ = 0, and was 
therefore omitted here. Summing over the channels ex- 
plicitely and expanding the distribution functions to lin- 
ear order in ST, we find 



ttGq 



^^ = ^TG , (5.23) 
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in agreement with the Wiedemann-Franz law, and 
A f dui f de 



An e i 
where 



4r / 2iri 



2n 



l-(Lu,e)eK H (u), (5.24a) 



i/S/l^>"*»- < 5 - 24b » 



A 



h L g 2 R + hRgj 

(gL + gfl) 3 



(5.25) 



when 



,16 



3 



(5.32) 



Equation (|5.23|) together with Eqs. 15.29|) gives the 
singlet channel part of the result for the thermal con- 
ductance reported in Sec. Ill Bl Evaluation of the triplet 
channel contribution is straightforward (as explained 
at the end of the previous subsection) and leads to 
Eq. l|2.22c|) . 



is the same form factor that appears in equation Ij5.13|l . 
and Z± are the following combinations of distribution 
functions: 



2-±(w, e) = ± [/ (e)A/(e - u) ± f (e - w)A/(e)] 



Note that at linear order 



A/(£) " 2T^7' 



(5.26) 



(5.27) 



where /o is the equilibrium Green function given in 
Eq. H3.8JI with temperature T and /i = 0. 

In Eqs. H5.24|) . we perform the integration over e: 



+ 



8Ttt 2 5 
7T 3 duj 



u> coth — ) 
27V 



4f- 2 i^w (5 - 28b) 



After substitution in Eq. I|5.24|l . the last term in the right 
hand side of Eq. (|5.28a|l gives a term proportional (after 
intergrating uj by parts) to the integral Eq. I|5.16[) appear- 
ing in the calculation of the electrical conductance, and so 
contributes to A« e i a term that obeys the Wiedemann- 
Franz law; the first term, on the other hand, gives a 
contribution proportional to the inelastic correction: 



An el 
A/Cj n 



TAG C 



A 



3e 2 

A Jin i 



(5.29a) 
(5.29b) 



where 

J' in. — 



1 f du> , d , u> r ^r>, . „„, 

In this expression, the final integration over u) can be 
done similar to Eq. I|5.16[) by closing the contour in the 
upper half plane, and the result can be expressed in terms 
of polygamma functions: 



J; , 



TFT 

1~2~ 



5i (2tttT) - gi 



2-ktT 



(5.31) 



C. Electrical conductance: Cooper channel 

The evaluation of the Cooper channel correction to the 
electrical conductance is based on the current formula 
()4.93|l . Writing the DC current as in Eq. Ij5.4|l and sum- 
ming over n, the contribution of the first term in square 
brackets in Eq. I|4.93|l gives the classical conductance Go, 
and the contribution of the other terms gives the Cooper 
channel correction which can be written as: 



AG Cooper — AG el + AG MT , 



(5.33) 



where 



AG e i = —A [J DoSl + J DoS2 ] , (5.34) 
AGmt = -BJ MT , (5.35) 
with the form factor A defined in Eqs. (|5.25(l and 



B = __gLgR_ _ A 
(gL + gfl) 2 



(5.36) 



which are characteristic of elastic and inelastic contribu- 
tions respectively, and we introduced 



rSdVj 2ttJ 2tt ^ 1 m| V ' 



m=±l 



A/ (| + e)j(^-e,u) \T>*\\ U ), 



2 7r d f du> f de v-^ o , . 



(5.37) 



m=±l 



(5.38) 



xA/(| + ej^|-e,o,)|2?£| 2 M, 



2nd f de f duj ( oj \ 

J2 Ci(2e)V*(u) 

m=±l 



x 9 ( 2 + e ) Im 



(5.39) 



The kernel J(e, oj) and the function A/(e) are defined in 
Eqs. (|4.91|l and l|5.2|l respectively. 
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To proceed further, we need the explicit form of the 
polarization operators in the presence of the bias voltage, 
as they determine the the kernel J(e : uj). In this case to 
zeroth order in interaction we have [see Eq. (|5.1U|) ]: 



Q = ±l 

where g±i are defined as 

_ gR 

g+i 



aV/2 



2T 



(5.40) 



(5.41) 



El + gR gL + gR 

Using the formula (obtained similarly to Eq. I|4. 101(1 ~) 

J de ReC m (2e) tanh ( 6 \ = Im^ (0) Q+* m ) , 

(5.42) 

with z m given in Eq. I|4.100[l . and the definition 

1 (\ — + i(u - mE% — aV)\ 

5 (2 + *i -)• 

(5.43) 

with ip^ the i-th polygamma function, from Eq. (|4.97|) 
we obtain 



2m - (q + (3)V 
AT 



m.a. f \ 

Xo M • 



illf (w) = ^2 g a g/3 coth 

(5.44) 

Here and below, the indices to, a and f3 are summed over 
±1; the quantities r* and E* z are defined in Eqs. 
and lHHIJ. Using again Eqs. l|512|) - l|513|) . Eqs. 
give: 



ImngH = -i^g QX ™' Q M, 



(5.45) 



Substituting the above results into Eq. H4.91[l we find: 

(5.46) 



m,a,0 



coth 



/ 2w- (a + f3)V S 



4T 



tanh 



2e-(3V 
AT 



What remains to be done, are the two integrals over to 
and e in Eqs. (|5.37|) - (|5.39l) : the latter can be evaluated 
exactly, while the former only approximately. Due to 
their similarity, we consider Jmt and JdoS% together in 
the next subsection, deferring the calculation of JdoS 2 to 
a later subsection. 



1. Maki-Thompson and D0S1 parts 

Substituting Eq. 1(5. 46|) into Eqs. 1(5. 37|) and using 
Eq. 1(5.42(1 together with the identity: 



-\C m \ 2 (2e) =ReC m (2e), 



(5.47) 



the result of the integration over e is: 

m,n a, (3 



AdV J 2?r 



(5.48) 



coth **-{a + (J)V _ cQth 2»-V-a)V 



AT 



AT 



•Xa M [Xo M+Xo M ) 



where g a and x™' Q (u;) were defined in Eqs. 1(5.41(1 and 
1(5.43(1 respectively, and the subscripts to, n, a and /3 are 
summed over ±1. Moreover since 

-ReC(2e) = — ^- T -RcC m (2e), (5.49) 
o(r* ) 

the expression for JdoSx is found by replacing the Xo' s 
with Xi/(27tt»T) in Eq. I|5.48|l . For convenience, in both 
Jmt and JdoSi > we separate two contributions with dif- 
ferent dependences on gi._R (and we drop a contribution 
odd in lj which vanishes upon integration) : 



_ 2g £ g fl . 

J CI / \ n *J a 



gl+gR T 2 



gL + gR) 2 ° (gL + gR 



(5.50) 



where a = MT, Do Si. For a = MT we have 

lay y 2^ |Pa| h ^ x ° {u)x ° (w) 
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xa 



coth 



1 - coth I — I 



2T 



\2T J 



(5.51) 



J MT — 



J m,a n,0 



xa 



coth 



uj — aV 
2T 



coth (—) 
\2TJ 



(5.52) 



and, as explained above, the formulas for a = D0S1 are 
found by replacing Xo by Xi/{^ nT *T). 

Up to now, no approximation has been made. However 
the validity of the RMT for the metallic dots requires all 
the energy scales [T,V,1/t*,E%] to be much smaller than 
T c ; we can then neglect the energy dependence of the 
propagator and approximate 



e 



(5.53) 



with e defined in Eq. ^2.14[l . The last approximation 
is valid with logarithmic accuracy and can be verified by 
noticing that the dependence of the polarization operator 
on the relevant quantities is the same as in Eq. (|5.45(l 
(dropping "Im" on both sides of the equation) and using 
Eq. H4.105f) . Even with this approximation, the integrals 
can not be calculated in closed form and we must resort 
to further approximations valid in different limits. In 
particular, we will consider the low and high temperature 
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regimes - the transition between the two occurs at T ~ 
1/t*, with t* defined in Eq. I|4.83f) . 

In the low temperature limit T <C 1/r*, we can use in 
Eq. (|5.43|) the following asymptotics for the polygamma 
functions: 



ImV (0) [\ 



Im^ (1) 



4nT 



4ttT 



arctanr„a; 



-AttTt * 



1 + (r*w) 2 



(5.54) 
(5.55) 



We further distinguish two cases: "low" and "high" volt- 
age, when V <C and V ^> T respectively; it is evi- 
dent that the two conditions are not mutually exclusive 
and the results derived below must agree with each other 
at intermediate voltages. 

If the voltage is small, V <C 1/r*, we can further ex- 
pand x"' Q (w) m T *y and t„cj as well (as \io\ is itself lim- 
ited by \V\) to obtain 



n,a, \ 2 t*(u - aV) 

n \ zj / 

Moreover in this approximation we have 



(5.56) 



so that in this regime 



, 1-(t,^) 2 



1 + (r*£|) 



(5.58) 



and we only need to calculate Jmt- 

Inserting Eq. I|5.56[) into Eq. (|5.51|) and performing the 
integral over u> we find 



71 



6(1 + (r*E* z )^ 



(5.59) 



where we used the identity 
8 



dy 



dx [coth(a; — y) — coth(x + y)] (x — y ) 



2 ^TT 2 
8y + ~- 



(5.60) 



By repeating the above steps for Eq. I|5.52[l we find 



J 



AIT 



e(l + (t*E* )2) 



2 r 



(5.61) 



Using the results l|5.59|l and l|5.61|l in Eq. (|5.5U|I and then 
Eq. H5.35f> . we arrive at Eqs. I|2.15|l . 

At high voltages V T (but still low temperatures 
T <C 1/t*), we can approximate the hyperbolic cotan- 
gents in Eqs. I|5.51|l - H5.52|) with their zero-temperature 



limit: coth ^ — ► sgnw. Performing the differentiation 

with respect to V and a partial integration, and using 
the approximation in Eq. (|5.54(l , we obtain 



J MT — 



- y 



n,m — ±1 



v 



t* arctan r* (uj + V — nE* z ) 
l + r|(w -V -mE* z ) 2 



(5.62) 



The result of the w-integral can be expressed in terms of 
dilogarithms: 



E h(r*V,T*V + mr*E* z ), 



(5.63) 



J 2 



m— ±1 



T MT 



JM 

MT ~ 2 1 As 2 



arctan t*(U — mE* z ) 



\m=±l 



The second line is found in the approximation (|5.54|l by 
direct comparison of the definitions (|5.51[1 - H5.52[1 and the 
function h is defined as: 

h ( x , y) = {f( x , y) - f(-x, y) + f(y, x) - f(-y, x)}, 

(5.64) 



with 



f{x,y) = Li 2 



x + y + i 
2(2/-*) 



Li 2 



x + y + i 



+ In (1 + ix + iy) 



In- 



2(2/-*) 



hi 



22/ 

y — x + i 



2>J 



(5.65) 



The approximate expression l|2.17|l for h is found as fol- 
lows: we first notice that in Eq. I|5.62l) we can restrict 
the integral to the region (0,V) (if we multiply by 2); 
then we shift the integration variable: u — > u> + V + mE* z 
and we finally neglect the dependence of the numerator 
on u>. The resulting integral is straightforward and gives 
the formula l(2~T7Jl . 

We now turn to the high temperature regime T 3> — . 
In this case, we can neglect 1 / (Attt^T) in the argument of 
the digamma function in the definition (|5.43(l ; then using 
the identity 



-lm-0 (o) I =tanhw, 

7T \ 2 7T / 



(Imw = 0), (5.66) 



we have 



■ tanh — 



Xo M 



aV — mE% 



AT 



(5.67) 



Substituting the above approximate expression into 
Eqs. (|5.51|) - (|5.52|) and performing both the differentia- 
tion with respect to V and the integration over oj we 
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obtain 



so that 



t1 

'J AIT 



1 7T J 1 /£* 

— < 1 — — coth —Z- 

2 e 2 I 2 V 2T 



(5.68) 



ci 



2T 



+ ci 



£7? 



u 



2T 




where 



c„(x) 



(a; cothx) . 



(5.69) 



(5.70) 



Since JdoSx is smaller than Jmt by a factor of (t^T)" 1 
[see Eqs. (|5.49|) ] , it can always be neglected in this regime. 

This concludes the calculation of the Maki-Thompson 
and "ReC 2 " contribution and hence the derivation of the 
results for AGmt reported in Sec. [HI In the next subsec- 
tion we present the calculation of the remaining "DoS" 
correction for completeness. 



D0S2 part 



In this subsection we evaluate Jdos 2 > Eq. 
the relevant limits. Using Eqs. (|5.2I) . I|5.40|l and I4.86|l . 
integration over e can be done similarly to Eq. (|4.101|l by 
closing the integration contour in the upper half plane. 
The result is 



1 ^-^ d fda; 
2 = 16t^T E ^ ga qvJ2tt C 



m,o>,{3 



2co-(a + (3)V 
4T 



x Re{[r'^(-^)-r' Q H]^M}, (5.71) 



where 



afr ^_ 1„, (1) (l , £ + t(w-m£7J -aVY 



4ttT 



(5.72) 

and the subscripts m, a and (3 are summed over ±1. In 
Eq. H5.71J1 , the approximation H5.53(l amounts to neglect- 
ing the imaginary part of the propagator in comparison 
with its real part, so that 

1 ^-^ d fdu ( Id \ 

Jd ° s ^i6^t E W ^ C0th l2rJ 

m— ±1 

xRe [r' + H - r ,_ H] Re^f(0). 

(5.73) 

In the low temperature regime T <C — , we can use the 
following asymptotic expansion: 



z>l 



(5.74) 



J^^Re^(0)Re£/^coth(^) 



dV 



i + i(u - aV - mE\ 



(5.75) 



The result of the integral can be expressed in terms of 
t/jW and using again the expansion (|5.74() we find 



Jd 



O.S-2 



1 E - 

2e J-i, 1 



(5.76) 



In the high temperature limit T ^> — , we can use 
Eq. I|5.67|l to rewrite JdoS 2 as 



ReX>?(0) v-^ d fdCj 
JdoS 2 « ^_ A V V /^coth2 W 
— oV 



16r*T 



m = ±l 



(5.77) 



sech 2 (cl> — V — mEz) ~ sech 2 (u) + V" — mEz) 



and integrating over w we obtain 

ir/e x (V-n,E, 



JD0S2 



16t*T 



E 

m=±l 



f'2 



2T 



(5.78) 



where C2 was defined in Eq. (|5.7U|I . We see that for 
t»T ~ 1, we have Jp s 2 ~ eJmt > |Jjmt|- However 
the prefactor A in Eq. (|5.34|) vanishes for reflectionless 
contacts and for this reason we only report the MT con- 
tribution in Sec. UT1 



VI. CONCLUSIONS 

In this paper we presented a quantum kinetic equation 
description of transport in open quantum dots with the 
inclusion of all the "universal" - in the Random Matrix 
Theory sense - interaction effects [sec Eq. (|1.5I) ]. While 
the effect of the charging energy was taken into account 
before^Sii the interaction corrections in the triplet and 
Cooper channels are considered here for the first time. 

The main result of the present work is that the triplet 
channel interaction can significantly affect the differen- 
tial conductance of the quantum dot, see Fig. [21 and in 
contrast with the singlet channel contribution, it is sen- 
sitive to the magnetic field, cf. Eqs. I|2.5|l and l|2.6|) - see 
also Fig-El F° r transparent contacts both the singlet and 
triplet channel corrections to the electrical conductance 
vanish; the non- vanishing Cooper channel contribution is 
unfortunately expected to be a negligible one in metallic 
dots, as discussed after Eq. (|2.14|l . 

In addition to the electrical conductance, we were able 
to calculate the thermal conductance by applying the lo- 
cal kinetic equation approach developed in Ref. Il6l For 
the thermal conductance we find that the Wiedemann- 
Franz law is violated by the interaction corrections [see 
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Eq. 1)2.20)1 ] . and we investigated the effect of magnetic 
field on the Lorentz ratio for contacts of finite reflection. 
The charge and triplet channel corrections to the ther- 
mal conductance also vanish for reflectionless contacts, 
and the Wiedemann-Franz law is not violated by the 
electron-hole channels in this case. 



We can now calculate the correction Scy to the specific 
heat and hence arrive at the expression for the energy 
density indeed 



8c\ 



l d 2 SQ du b 
dT 2 ~ W 



(B3) 
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Using the first of the above relations and after an inte- 
gration by parts we find 



Scy — 



d_ 

dT 



<1ujujNp(lo) b p (u))~b 9 {uj) , (B4) 



APPENDIX A: VANISHING OF St 
STEADY STATE 



IN THE 



We want to show that the properties ()4.43)1 and (|4.5U|) 
hold in the steady state. To this end, we notice that 
the second term on the right hand side of the definition 
(|4.22b|) always satisfies these properties, as one can verify 
using the same arguments that prove the corresponding 
properties (|4.42l) and 1)4. 46b)) of the collision integral St+ . 
Next we want to show that in the steady state 



[SQ] K_ = 0, 



(Al) 



after the average over the fluctuating field but before any 
limit and/or time derivative is taken. Indeed in 5Q only 
the contributions proportional to 8g^_ should be kept, 
since the ones proportional to Sg^ vanish after the aver- 
age. We note that the right hand side of Eq. (|4.12a|l can 
be rewritten as 



i[K + ;d t g K ] , 



(A2) 



by using the zeroth order (in the interaction) part of 
Eq. H4.19fl for g K , i.e. neglecting the collision integrals 
Sti and St2 - this is sufficient in the one-loop approxi- 
mation. In the steady state, dtg = and hence 5g+ = 0. 
This proves Eq. i)Al)l and therefore the properties 
and g30J). 



where we introduced the bosonic density of states: 

b a {u) = -Im^ln£ Q . (B5) 

7T 

Comparing Eq. I)B4)1 to the second relation in Eq. i)B3)) 

and generalizing this result to an arbitrary distribution 
function we have 



Ub = / dio uj 
la 



N p (u>)b p (uj)-N g (uj)b 9 (uj) . (B6) 



It is straightforward to prove that this coincides with 
Eq. H4.49[> [by substituing in the latter Eq. (|4.39|) and 
noticing the different limits of integration] . 




+ 




APPENDIX B: THERMODYNAMICS 

The effects of interaction on the thermodynamics can 
be obtained by calculating the leading contribution to the 
themodynamic potential f2; this is given by the sum of 
the so-called ring diagrams of Fig. |S][see e.g. Ref.||. This 
sum can be calculated in the Matsubara representation 
and after standard analytic contnuation the correction 
Sil is found to be: 



(50 = / coth 

2vr 2 



— I Im In 

IT) 



1-2£JI£H . (Bl) 



Using expression H4.29J1 for the polarization operator and 
the definitions 1)4. 7JI and 1)4.32(1 . we can rewrite this as: 



f duj 






\nC 9 




- coth 


— ) Im 


J 2^ 


2 


k2TJ 





(B2) 
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FIG. 8: Leading singular contribution to the thermody- 
namic potential. The shaded box corresponds to 2E C , defined 
through the two particle vertex [see e.g. Ref . l24| : the solid 
lines are coherent parts of the electron Green functions. 
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